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Abstract Gamma-Ray Bursts (GRBs) are the strongest explosions in the Universe, 
which due to their extreme character likely involve some of the strongest magnetic 
fields in nature. This review discusses the possible roles of magnetic fields in GRBs, 
from their central engines, through the launching, acceleration and collimation of 
their ultra-relativistic jets, to the dissipation and particle acceleration that power their 
y-ray emission, and the powerful blast wave they drive into the surrounding medium 
that generates their long-lived afterglow emission. An emphasis is put on particular 
areas in which there have been interesting developments in recent years. 
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1 Introduction 


Gamma-ray Bursts (GRBs) are among the most extreme objects in the Universe. 
They are the most luminous cosmic explosions, and therefore serve as beacons at 
the edge of the visible Universe that can be used as cosmic probes. GRBs provide 
short timescale insight into end-stage stellar evolution, and serve as probes of extreme 
physics such as strong gravity, very large densities and magnetic fields, extremely 
energetic particles, and relativistic bulk motions. They are also promising sources of 
high-energy neutrinos and gravitational waves. 

GRBs can be roughly divided into two main sub-classes: (i) Long-duration (> 2 s) 
soft-spectrum bursts that are found in star-forming regions and are associated with 
broad-lined Type Ic supernovae, implying a massive star progenitor, which is most 
likely low-metalicity and rapidly rotating near this cataclysmic end of its life, and 
lives in a gas-rich environment not far from where it was formed. In order to pro¬ 
duce a GRB, the central engine must drive a strong relativistic jet that bores its way 
through the stellar envelope and produces the GRB well outside of the progenitor 
star; (ii) Short-duration (< 2 s) hard-spectrum bursts that are thought to arise from 
the merger of a binary neutron star system (or a neutron star and a stellar-mass black 
hole) that emits gravitational waves as it inspirals and coalesces, producing a cen¬ 
tral engine driven jet. Such systems live in low density environments, possibly with 
a prior supernova kick that pushed them into the outskirts of their host galaxies. A 


third sub-class, whose importance was realized only recently (Soderberg et al. 2006 

Campana et al. 2006 

Liang et al.|2007; [Virgili et al.||2009| Bromberg et al. 2011 

2012 Nakar and Sari 

2012), involves low-luminosity GRBs, whose overall isotropic 


j y,iso 

single-peaked lightcurve, and their vF v spectrum typically peaks at a lower than av¬ 
erage photon energy (usually E v < 100 keV). While observed rarely, because of their 
low luminosity, they are the most numerous group in nature (in terms of their rate per 
unit volume). They most likely do not arise from the same emission mechanism as 
regular long GRBs (e.g., |Bromberg et al.|2012||Nakar and Sari|2012| ). 

The phenomenology of GRBs is generally separated into two observational phases: 
the prompt emission and the afterglow. These two phases are traditionally differenti¬ 
ated largely based upon instrumental measurement methods, but they do seem to also 
be physically distinct - they arise from different emission mechanisms and occur at 
different distances from the central source. However, the dividing lines between the 
prompt emission and the afterglow have blurred in recent years. In the standard Fire¬ 
ball model (e.g., FiranT999, 2004| Kumar and Zhang 2015), the prompt emission 
(i.e. the burst of y-rays) is due to dissipative internal shocks within the outflow, while 
the long-lived broadband afterglow is the result of the jet driving a strong relativistic 
forward shock into the surrounding medium as it decelerates and transfers its energy 
to the shocked external medium. Unless the outflow is highly magnetized when it 
is decelerated by the external medium, this deceleration can occur through a strong 
reverse shock that results in a bright optical flash, which is also sometimes detected 


(e.g., Akerlofetak|1999 JSari and Piran|1999| Mundell et ah 2007b| Racusin et al. 


|2008[ Vestrand et al.||2014 1 just after the onset of the prompt emission, and decays 

largely independently of the forward shock emission. 
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The field of GRBs is relatively young, with several revolutions in our understand¬ 
ing of these objects thanks to new observations over the last two decades. The role 
of magnetic fields in GRBs is relevant to many topics in this field. They affect the 
properties of the compact object (neutron star and/or black hole) that powers the cen¬ 
tral engine, and how it launches the jet. Magnetic fields may also play an important 
role in the acceleration and collimation of the relativistic jets in GRBs, as well as in 
their composition. They can contribute to the energy dissipation and particle acceler¬ 
ation that powers the prompt GRB emission, and may play a key role in its emission 
mechanism. A strong magnetic field can suppress the reverse shock and its emission. 
During the afterglow, the amplification of the weak magnetic field in the external 
medium by the afterglow shock and its subsequent behavior in the shocked external 
medium downstream of the shock play a key role in the particle acceleration by the 
shock and in the shaping of the afterglow emission. In fact, it seems hard to find any 
important part of GRB physics where magnetic fields might be safely ignored. 

In this review, we explore the evidence for extreme magnetic fields in GRBs, and 
how magnetic fields are intertwined with our understanding of the mechanisms that 
produce the relativistic jets that power these objects. As we cannot cover here all 
of the relevant topics in detail, we have instead chosen to focus on specific topics 
in which there has been recent progress (see e.g. Piran|2005j for an earlier review). 
First, in <)2]a brief overview is given on the progenitors of both long and short GRBs, 
with the main thrust being devoted to the possible role of millisecond magnetars - 
newly born, very rapidly rotating and highly magnetized neutron stars - a topic that 
has recently received a lot of attention in the literature and in the GRB community. 
Next, £[3] discusses the dynamics of GRB jets. It starts with long GRB jets as they 
bore their way out of their massive star progenitors, and then moves on to discuss 
more generally the possible role of magnetic fields in the acceleration and collimation 
of GRB jets, both in steady-state and highly time-variable outflows, as well as in 
the interaction of the jet with the external medium and the reverse shock. In <[4] we 
discuss the role of magnetic fields in the dissipation and radiation that power the 
prompt y-ray emission, and what GRB observations can tell us about the conditions 
within the emitting region. Finally, |j5]is devoted to the role of magnetic fields in the 
afterglow. It focuses on their effects on the afterglow and reverse shock emission and 
their polarization, and how this can teach us about the magnetic field structure in the 
GRB outflow and its amplification in the afterglow shock as well as its structure and 
possible decay further downstream of this shock. Our conclusions are discussed in 


2 GRB Progenitors, Central Engine, and the Role of Magnetars 


Long duration GRBs are associated with Type Ic supernovae (e.g., Woosley and 
Bloom 2006) ), which directly relate them to the death of massive stars stripped of their 
hydrogen and helium. This supports the popular Collapsar model (Woosley 1993 


MacFadyen and Woosley 1999j ) in which a central engine lunches a relativistic jet 
that penetrates the stellar envelope and powers the GRB. Typically, within the Col¬ 
lapsar models the central engine is considered to be a an accreting newly-formed 
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stellar-mass black hole at the center of the progenitor star. The most popular model 
for short duration GRBs features the merger of two neutron stars in a tight binary 
system (Eichler et al. 1989j >, which may again form a black hole surrounded by an 
accretion disk as they coalesce. Therefore, in both long and short duration GRBs, 
despite their very different progenitors, the central engine that is formed during the 
explosion and launches the relativistic jets might still be similar in nature - accretion 
onto a newly formed black hole. 

An attractive alternative possibility that has gained popularity in recent years is 
that GRB central engines may involve magnetars ( Usov|1992j Duncan and Thomp 


son 1992 Bucciantini et al. 2008 Dall’Osso et al. 20091 - highly magnetized neutron 
stars with surface magnetic fields of order B ~ 10 13 G, which in this case are newly 
born and very rapidly rotating, with ~ 1 ms periods (and hence dubbed millisec¬ 
ond magnetars). In this model the main energy source is the neutron star’s rotational 
energy, and a very strong magnetic field is needed for a rapid extraction of this ro¬ 
tational energy and to channel it into a relativistic outflow. A rapidly rotating neu¬ 
tron star may arise in the collapse of a rotating stellar core and the magnetic field 
can be amplified in this collapse (e.g., Duncan and Thompson [1992| ). As the magne- 
tar’s energy is naturally extracted in the form of a Poynting flux (though this flux is 
initially not significantly collimated) it naturally leads to a magnetically dominated 
outflow. Collimation of the outflow into a narrow jet may, however, be facilitated by 
the interaction of the outgoing strong magnetohydrodynamic (MHD) wind with the 
progenitor star’s envelope (e.g.,|Bucciantini et al.|2007||2008 2009;|Bromberg et al. 
|20T4| ). 


After the first ~ 10 — 100 s or so from the neutron star formation, the neutrino- 
driven winds subside and the baryon loading on the MHD wind significantly de¬ 
creases. As a result the initial wind magnetization parameter Oq significantly in¬ 
creases and becomes 3> 1 (the magnetization parameter a is the Poynting-to-matter 
energy flux ratio, or proper enthalpy density ratio). At this stage, the neutron star 
spin-down and its associate luminosity are approximately given by the magnetic 
dipole in vacuum formula (which also approximately holds in the force-free regime; 
Spitkovsky|2006| l, 


B 2 R 6 Q$/(6c 3 ) 


= u » ^oMUi ) _ Eo 

[l + 2B 2 R 6 Qyt / (6Ic 3 )] 2 t 0 (l+t/t 0 ) 2 


^L 0 x 


(t < to) 


m 


where I is the neutron star’s moment of inertia, R is its radius, B is the surface dipole 
magnetic field at the pole, Qo is the initial angular velocity, Eo ~ 4 / £2/ is the initial 
rotational energy, Lq = Eo /to is the initial spin-down luminosity and 

3/c 3 , ( B \ “ 2 / P \ 2 

'» = ^= 2x10 (tcPg) (is;) s (2) 

is the initial spin-down time (using typical values of R ss 10 km, I « 10 45 g cm 2 ). 
This spindown luminosity initially (at t < to) has a plateau at Lq, and then (at t>to) 
falls off as t 2 . Both Lq and to can be tuned with the proper choice of the initial 
angular velocity £2o and the magnetic field B. With a choice of B ~ 10 15 ' 5 16 G one 
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can arrange to fit the prompt duration in which case the magnetar is invoked to power 
the prompt GRB. With a lower magnetic field of order 10 14 5 ~ 15 G, to is of order 
several thousand seconds, comparable to the duration of the plateau phase in some 
X-ray afterglows. 


The magnetar model gained a lot of popularity with the discovery by Swift of 
plateaus in the X-ray afterglow light curves of many GRBs (Nousek et al. 2006[ 
Zhang et al. 2006[>, whose shape resembles the overall shape of Magnetar’s spindown 


luminosity (Tr oja et al.||2007| Dall’Osso et al. 2009 201 1| Rowhnson et al.||2013 


20141. Somewhat surprisingly, even though this tentative evidenc^jfor magnetar-like 


activity was obtained for the afterglow phase, it was interpreted in the community 
as evidence for a magnetar operating as the main source of energy for the prompt 
emission as well. Both interpretations face some difficulties. 


The magnetic field needed to produce the prompt emission is larger by about one 
order of magnitude than the one observed even in the strongest magnetars. This may 
not be that puzzling as there is ample evidence of magnetic field decay in magne¬ 
tars (Dall’Osso et al. 2012} , and the observed magnetars are typically a few thou¬ 
sand years old. It is possible that the magnetic fields of newborn magnetars are large 
enough. 


A more serious problem concerns the energy budget. The rotational energy of 
a typical neutron star, even when rotating at breakup velocity, is at best marginally 
sufficient to power the most powerful GRBs ( |Cenko et al.|2010} . This is especially so 
if we also take into account the efficiency of converting this rotational energy into the 
prompt flux of y-rays. Of course, magnetars could still power less energetic GRBs. 
However, this would require one to invoke two kinds of central engines, as a different 
energy source would be needed to power the most energetic GRBs. 


Even more perplexing is the situation concerning the longer duration plateaus 
in the afterglow light curves. Here, the needed values of magnetic fields are indeed 
typical for those arising in the observed magnetars, and the overall energy budget is 
reasonable as well. However, another question arises: if a low magnetic field mag¬ 
netar has powered the afterglow plateau, then what has powered the prompt GRB? 
Can a magnetar fire twice? The simple answer, according to Eq. ([lj, is no. This is as 
long as the magnetic field remains constant during the slowdown time scale. How¬ 
ever, one can come up with a fine-tuned model in which the magnetic field decays 
on a timescale shorter than tg. In this case the duration of the magnetar activity is 
not determined by to, but by the magnetic field decay time. Once the magnetic field 
has decayed, a second slower magnetar phase appears with a new to- Overall such 
a model requires a the magnetic field that is extremely large initially, leading to the 
prompt emission and then it decreases, just at the right time (and before all the ro¬ 
tational energy is exhausted) to a lower level in which the weaker magnetar powers 
(using the remaining rotational energy) the afterglow plateau. 


1 These plateaus have several alternative explanations, which are at least as compelling as the magnetar 
explanation, such as promptly ejected slow material that gradually catches up with the afterglow shock 


i Nousek et al. 2006 (Granot and Kumar 2006), time-varying afterglow shock microphysical parameters 

Granot et al. 2006 

, viewing angle effects jEichler and Granot 2006 ), or a two-component jet l|Peng et al. 

2005 Granot et al. 

20061. 
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An alternative option is as follows: in the first ~ 10 — 100 s or so after the for¬ 
mation of the neutron star, the strong neutrino-driven winds cause a large baryon 
loading on the MHD wind that prevents the formation of a very high initial magne¬ 
tization ((To < 10 100). Therefore, during at least part of this time the spindown 

luminosity can significantly deviate from the form Eq. 0 - the formula for a mag¬ 
netic dipole in vacuum - and may in fact be significantly higher, and closer to the 
result for a magnetic monopole in vacuum, since most of the magnetic field lines are 
opened by the strong baryon loading (e.g.,[Metzger et al. 2011) >. This can increase the 
spindown luminosity by a factor of ~ ( R // R) 1 ~ 10 1 * '^7where R/ =c/£2q is the initial 
value of the light-cylinder radius. However, as in the early magnetic field decay sce¬ 
nario mentioned above, also this solution would require fine tuning in order to extract 
just the right amount of rotational energy over just the right timescale. Moreover, in 
this case the bulk of the large amount of energy that is released on the timescale of 
the prompt GRB is given to relatively low-Oo baryon-rich material, which could not 
attain sufficiently large asymptotic Lorentz factors that are needed to power a GRB. 


Another possible solution to this problem was suggested recently (Rezzolla and 
Kumar|2015| Ciolfi and Siegel|2015[> - the “time reversal model”, which postulates 


magnetar activity for the plateau but an accretion disk for the prompt phase. Ac¬ 
cording to this model, first a magnetar with to ~ 10 4 s is born and launches a fast 
MHD wind whose interaction with slower matter that was ejected earlier produces 
the afterglow plateau. In this scenario the magnetar is a supramassive neutron star, 
i.e. supported against gravitational collapse by its very fast rotation, so once it spins 
down significantly it collapses to a black hole, and an accretion disk that forms dur¬ 
ing this collapse powers the prompt emissionj^jln spite of this reversed time sequence 
the plateau is observed after the prompt emission because it involved the interaction 
of the winds and this phase introduces a time delay (in the observer frame). While 
it is appealing, this model requires the formation of a disk during the collapse of the 
supramassive neutron star. However, [Margaht et al. ( 2Q15[ > have recently argued that 
this is impossible. 

A different possible solution is if the prompt GRB is powered by the energy in 
a strong initial differential rotation ( Kluzniak and Ruderman|1998| ). The strong dif¬ 
ferential rotation winds-up strong toroidal magnetic field loops, which are buoyantly 
pushed out of the neutron star surface and power the prompt GRB. This lasts until 
they exhaust all of the differential rotation energy on the timescale of the prompt 
GRB emission of long duration GRBs. The rotational energy of the remaining uni¬ 
formly rotating neutron star could then power the plateaus on its longer magnetic 
dipole spindown time to. This might possibly work for long duration GRBs that are 
not too energetic (as the energy in differential rotation is somewhat lower than the 
total rotational energy). 

Recently, it was suggested that millisecond magnetars might also be at work in 
short duration GRBs (e.g., FanmidXu]|2006 Rosswog||2007 Metzgeretak|2008 


IRowlinsorP et al. 20T0j |Bucciantini et al. 20121 Rowlinson et al. 20131 |Fan et al.| 


2013a[ Gompertz et al.|2014ji. Newly formed millisecond magnetars were suggested 


2 This scenario is rather similar to the “supranova model” that was suggested much earlier Vietri and 

|Stella|199B|[T999). 
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to produce the extended emission seen on a timescale of ~ 10 2 s in some short GRBs 
(e.g. Metzger et al. 2008 Buccianti ni et al.|2012| l. In this picture the initial short GRB 
spike may arise from short-lived accretion following the merger. The extended emis¬ 
sion is driven by the spindown power that is released over ~ 10 2 s, and takes several 
seconds to break out of the surrounding mildly relativistic material that is ejected 
quasi-isotropically during the merger. In this scenario, however, the jets that power 
the short GRB itself (the initial hard spike) are launched within the first second or so 
after the formation of a newly-born millisecond magnetar, when the neutrino-driven 
wind is very vigorous, the magnetization is low (<7o < 1), and the baryon loading 
is very significant near the star and the inner accretion disc where the jet may be 
launched. It is therefore unclear whether in this case the jets could eventually reach 
a high enough Lorentz factor to produce the GRB. Moreover, in a binary merger the 
neutron star is formed extremely rapidly rotating (near breakup), and its rotational en¬ 
ergy of a few x 10 52 erg is eventually injected into the afterglow shock. This should 
naturally produce a bright afterglow emission while the observed afterglows of short 
GRBs (either with or without an extended emission) are typically much dimmer than 
those of long GRBs. While this might in part be attributed to a lower external den¬ 
sity on average, this cannot fully account for the dimmer afterglows of short GRBs 
essentially over the entire broad-band spectrum, from radio to GeV energies. 


3 Jet Propagation and Dynamics 

The question of the jet composition is still a major issue in our understanding of 
GRBs. It affects the location of the emission site, the mechanism of the emission 
and the particle acceleration. There are two main possibilities that are commonly dis¬ 
cussed in the literature: a hydrodynamic jet and a Poynting flux dominated jet (for the 
jet dynamics discussed here we do not make the distinction between baryonic and e ± 
pairs particle content). The main advantage of a hydrodynamic jet is fast and robust 
acceleration, which allows the jet to reach very high Lorentz factors relatively close 
to the central source. Magnetic acceleration, on the other hand, is slower and less 
robust. However, hydromagnetic jet launching implies dynamically strong magnetic 
fields near the central source, which can naturally avoid excessive baryonic loading 
into the central part of the jet, and thus allow it to reach large asymptotic Lorentz fac¬ 
tors far from the source. The required very small baryon loading is hard to naturally 
achieve in a purely hydrodynamic jet. Here we focus on the possible role of magnetic 
fields in the jet dynamics and propagation, while keeping in mind these two main 
options for the jet composition. 


3.1 Jet Propagation in the Stellar Envelope 

In order to produce a regular GRB, a collapsar jet needs to successfully break out 
of its progenitor star. After it breaks out, the jet can accelerate freely and eventually 
generate the observed y-ray photons far from the star in a region where they can 
escape (see jQ. Before it emerges from the stellar surface, the jet propagates inside 
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the star by pushing the stellar material in front of it, forming a bow shock ahead of the 
jet. The stellar material that crosses this shock is heated and forms a cocoon around 
the jet, which in turn applies pressure on the jet and collimates it. The collimated jet 
propagates at a different velocity than a freely expanding jet. It continuously injects 
energy into the cocoon through a slower moving head that forms at the front of the 
jet. The head dissipates the jet’s energy and channels it into the cocoon. Therefore, 
the continuous propagation of the jet through the star depends on the supply of fresh 
energy from the source. If the engine stops injecting energy, the head will essentially 
stop propagating once the information about the energy cutoff will reach it, and the 
jet will fail. The breakout time, t b , is defined as the time of the engine shutoff for 
which the information about the shutoff reaches the jet’s head when it is at the edge 
of the star. If the engine stops working at a time the head will “feel” this cutoff 

while it is inside the star and will stop propagating. In this case the jet will not break 
out and it will not produce a regular GRBrl Since the information travels outwards at 
very close to the speed of light, the breakout time is related to the time at which the 
jet’s head reaches the outer edge of the star through 

f _ f Rt dz R*_R* 1-ft, 

b J 0 P h (z)c c ~ c p h ' 

where j3 h (z)c is the instantaneous velocity of the jet’s head at a distance z from the 
central source (along its symmetry axis, using cylindrical coordinates), and /3 h c is its 
average velocity. 

Simple analytic solutions to Eq. ([3]) can be obtained in two limits: (i) a non- 
relativistic limit, characterized by a proper speed (in units of c) = I/j3 h <C 1 (where 
r h is the head’s Lorentz factor) in which t b ~ R, //3 h c, and (ii) the relativistic limit, 
characterized by wj, 1, in which f b ~ R.,/2ff:. The transition between these two 
limits occurs when t b ~ R t j c, which according to Eq. ([3j> corresponds to j3 h ~ 1/2. 
The jet’s head is initially sub-relativistic, but it accelerates in the steep density profile 
inside the star (q = —d log p /cl log r > 2). Therefore, if the jet becomes relativistic at 
some radius, R iel , where Wh ~ 1, then it will remain so until it will break out. 

In order to calculate the breakout time ones needs a proper model for the propa¬ 
gation of the jet in the star. Such a model needs to consider the evolution of the jet 
and the cocoon self-consistently, as they affect one another. The propagation velocity 
of the head is determined by its cross section, which is set by the collimation of the 
jet. The head’s velocity, in turn, controls the energy injection into the cocoon, which 
determines the collimating pressure. The dynamics of this system can be described 
in a relatively simple way in two extreme cases of a purely hydrodynamic jet and a 
purely electromagnetic (Poynting flux dominated) jet. 

3.1.1 The Breakout Time of Collapsar Jets 

Close to the injection point the jet’s internal pressure, pj, is much larger than the co¬ 
coon’s pressure, p c . Therefore, initially the jet material expands freely until the col- 

3 A failed jet produces, most likely, a low-luminosity GRB when a shock wave generated by the dissi¬ 
pated energy breaks out from the seller envelope. 
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limation point where the jet’s pressure equals the cocoon’s pressure, p i = p c . Above 
this point the jet is collimated by the cocoon’s pressure, and its behavior depends on 
its magnetization. 

In the hydrodynamic case, the collimation of the jet leads to the formation of a 
collimation shock at the base of the jet (e.g., Bromberg and Levinson |2007j >. Above 
this point the jet maintains a roughly cylindrical shape due to a relatively uniform 
pressure in the cocoon (e.g., |Zhang et al.|2003||Matzner|2003j|Bromberg et al.|201 l| l. 
The jet material remains relativistic with a roughly constant Lorentz factor F ~ 1 / Qq, 
where Oo is the jet injection (or initial) half-opening angle. At the head of the jet 
the relativistic jet material decelerates abruptly through a strong reverse shock. Since 
the jet is roughly cylindrical upstream of the reverse shock, the width of the head - 
its cylindrical radius rj and its corresponding cross-section Ej = nrj, are set by the 
width and cross-section of the jet at the collimation point, which are shown to be 


v- 2 Mn 
Ej = nr ~ 

1 4 cp c 


(4) 


where Lj is the (one sided) jet luminosity. The velocity of the head was shown 
( |Matzner|2003| |Bromberg et al.|201 l| l to follow 


Ph 


Pj 

1+L-V2 ’ 


(5) 


where the dimensionless parameter 

£_ Pi h i r f ^ L i 

Pa EjP a ( 3 


(6) 


represents the ratio between the energy density of the jet (Lj/Ejc) and the rest-mass 
energy density of the surrounding medium ( p a c 2 ) at the location of the head. Here 
hj = 1 + 4pj/pjC 2 is the dimensionless specific enthalpy of the jet material just up¬ 
stream of the termination shock at the base of its head. 

|Bromberg et al.| ( |201 lj ) have obtained approximate analytic expressions for the 
propagation velocity of a hydrodynamic jet, and demonstrated that for typical stel¬ 
lar and jet properties, the jets head propagates at a velocity that is at most mildly 
relativistic. Therefore, in this case the solution to Eq. ([3]» is in the transition region 
between the relativistic and the non-relativistic limits, and can be approximated fol¬ 
lowing [Bromberg et al.j ( |2015| ). In order to obtain a useful analytic solution, they 
approximated the exact integration (shown by the dashed blue line in Fig. U} by: 


L.hyd — 6.5 




(7) 


where L re i is the transition luminosity between a non-relativistic and a relativistic 
case: 


Trel — 1.6 X 10 $0.84 


3 — < 


|\ 7/5 /5-| 


0.5 J 


2.5 


4/5 


7~£ 

4.5 


15/2 


ergs 


(8) 
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Fig. 1 The breakout time, f b , as a function of L calculated for a jet with an opening angle 9j = 7°, and a 
star with a mass M. = 1 5M .. radius R. = 4R@ and a power-law density profile p °= r 2,5 . The gray solid 
curve tracks the exact integration of Eq. {7}. the red and magenta lines show the analytic approximation 
for the non relativistic and the relativistic cases respectively. The dashed blue line follows the smoothed 
analytic solution for f b from Eq. jTJ. (This figure is taken from Bromberg et al. 2015}. 


As canonical parameters we have used here a stellar mass of M* = 15M 0 , a stellar 
radius R, = 4R & and we assume a power-law density profile: p, r “= with ^ = 
2.5. Hereafter we measure masses and radii in units of solar mass and solar radius 
respectively and use the subscript V to denote properties of the progenitor star. For 
all other quantities we use the dimensionless form A x = A j 1 O' measured in c.g.s 
units. For a typical collapsar (one sided) jet luminosity of Lj ~2x 10 49 erg s 1 , and 
injection angle of 0j = 7° the corresponding breakout time is fb( L,ei) —9 s. 

In a Poynting flux dominated jet the situation is different. The cocoon’s pressure 
is typically strong enough to collimate the jet before it looses causal contact with 
the axis. In this case the poloidal magnetic field is comparable to the toroidal field 
in the comoving frame of the flow, and shocks are inhibited. This leads to a smooth 
transition of the jet material from a free expansion state, close to the engine, to a 
collimated state. The jet material remains in a strong causal contact also above the 
collimation point. Therefore as it approaches the head, it does not shock. Moreover, 
it can be shown (Bromberg et al. 2014| > that under these conditions the jet’s proper 
velocity uj is approximately equal to the ratio of its cylindrical radius r and the 
light cylinder radius r L : uj ~ r/r L , and the same also holds at the jet’s head, iq-, ~ 
r h /r L . Therefore, the jet material gradually decelerates and becomes narrower as it 
approaches the head until at the head its velocity matches that of the shocked stellar 
material just behind the front tip of the bow shock. This deceleration and narrowing 
of the jet towards its head is assisted by the fact that the cocoon’s pressure becomes 
larger closer to the head, as the bow shock is stronger there. This results in a jet head 
that is much narrower than in the hydrodynamic case and therefore leeds to a much 
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faster propagation speed, where the head’s proper speed Wh is given by ( Bromberg] 
let al.|2014| ): 

a 1 / 4 5 (s h « 1) , 


where the dimensionless quantity 


a = 



npa^rl 


(9) 


( 10 ) 


is the ratio of the jet’s magnetic pressure at the light cylinder and the ambient medium’s 
rest mass energy density near the head Q 

Therefore, a Poynting flux dominated jet becomes relativistic at a radius R r deep 
inside the star, even with a modest power (Bromberg et al. 2014[ ): 


R , 

-^~1.4x 10“ 
R * 


49 3M*, 15M e «* i 4r q r L1 


R 


3-^ 

0.5 


i/« 


( 11 ) 


This implies that here only the relativistic asymptotic solution « a 1 / 6 ) is relevant. 
The corresponding breakout time is ( [Bromberg et al.|2015| : 


^b, 


mag 


, M . ,2/3 (0S_\ 

'49.3 m *,l5M 0 r L j l 3 _ ^ ) 


~o.8 lJ1 3 m! /3 


2/3 


( 12 ) 


This time is much shorter than the breakout time of a hydrodynamic jet with a similar 
luminosity. 

3.1.2 Observational Evidence for the Jet Breakout Time 


After the jet emerges from the stellar envelope it dissipates its energy at a large dis¬ 
tance and produces the GRB. On average, the overall behavior of the prompt emission 
does not vary significantly during the burst (the second half of the prompt emission is 
rather similar to the first one). This suggests that the prompt emission arises at a more 
or less constant radius and not in a propagating single shell. A single shell would have 
expanded by a factor of ~ 10—100 during the duration of a burst and it is unlikely to 
maintain constant conditions as it emits the prompt y-ray emission over such a wide 
range of radii. This implies, in turn, that the GRB activity follows the central engine’s 
activity ( Sari and Piran|1997] ), and that the GRB lasts as long as the central engine 
is active. Therefore, within the Collapsar model, the observed GRB duration (usually 
denoted by 790 , which measures the time over which the central 90% of the prompt 
photon counts are detected) is the difference between the engine operation time, t e , 
and the breakout time, f b , namely 790 =t e —t b (not accounting for the cosmological 
time dilation here). 

The breakout time essentially serves as a threshold time: a regular GRB is formed 
Bromberg et al. (2012]> have shown that in such a case one would 


only if t e > t } 


v 


4 This analysis does not account for 3D effects that can slow down the head’s propagation speed 

(Bromberg & Tchekhovskoy 2015, in prep). 






























12 


Granot et al. 


expect a plateau in the duration distribution of GRBs, dNcRB / dTgo, at durations that 
are shorter than t b . The logic behind this is as follow. At the time when the jet’s head 
breaches the edge of the star, it is already disconnected from the engine and cannot 
transmit information backward to the engine. In other words, the engine cannot “tell” 
when the jet breaks out of the star and we do not expect that t e and t b will be strongly 
related to each other. In fact, for a given t b we expect to have a distribution of engine 
activity times, where some are shorter ( t e < t b ) and some are longer (t e > t h ) than 
t b . In this case the probability of observing a GRB with duration Tgg is equal to the 
probability that the engine will work for a time t e = Tgo + f b : Egrb^o) = P e {t e = 
7i)o +t b ). This probability has a simple description in two limits: 


/ j grb(79o) ; 


Pe(t h ) (Tg 0 </„)> 

P e (Tgo) (790 >0- 


(13) 


Now, if there is a dominant population of GRBs with a typical t b , then at short dura¬ 
tions /V;rb(71)o) —> P e (t b ) = const, we expect to get a plateau at durations Tgo <C t b . 

Figure[2 depicts the duration distribution, dNGRB/dTgo , of BATSI0(21OO GRBs), 
Fermi-GBN ^] (1310 GRBs) and .S - vv//0(8OO GRBs). To fit a plateau in each data set 
[Bromberg etal.| ( [2015] ) looked for the maximal number of bins that are consistent with 
a plateau at a confidence level < 95% (2crJ^] The extent of the best fitted plateaus is 
5-25 s in the BATSE data (7.19/4 X 2 /d.o.f ), 2.5-17 s in the Fermi -GBM data (10/5 
X 2 /d.o.f), and 1-20 s in the Swift data (15.85/9 x 2 /d-o.f). Bromberg et al. (2015) 


accounted for three free parameters in the fit: the height of the plateau and the two 
opposite ends of the plateau line. The differences between the maximal durations 
of the plateaus can be mostly attributed to the different sensitivity and triggering 
algorithms of the different detectors. 

At short durations, the plateau is concealed by the increasing number of non- 


Collapsar (“short”) GRBs having a typical duration of less than a few seconds (Bromb ;rg 
et al.|2013| . As non-Collapsars have on-average harder spectra than Collapsars (e.g. 


Kouveliotou et al7||1993|, the relative number of non-Collapsars can be reduced by 


choosing a hardness threshold (for each sample) and selecting only the events that 
are softer than this threshold. This should lead to a less prominent “bump” at short 
duration. If the plateau is indeed an intrinsic property of the (softer) Collapsars du¬ 
ration distribution, it should extend to shorter durations in a softer subsample. To 
examine this effect they selected in each sample all the events that are softer than the 
median hardness of long GRBs (Tgo > 20 s) in the sample (see |Bromberg et al.|2013| 
for further details). Figure [3] shows the duration distribution of the soft GRB sub¬ 
samples. The plateaus indeed extend to much shorter durations than in the complete 
samples, supporting our hypothesis. The extent of the best fitted plateaus is 0.4-25 s 


5 http://gammaray.msfc.nasa.gov/batse/grb/catalog/current/ from April 21. 1991 until August 17, 2000. 

6 http://heasarc.gsfc.nasa.gov/W3Browse/fermi/fermigbrst.html, from July 17, 2008 until February 14. 
2014. 

7 http://swift.gsfc.nasa.gov/archive/grb_table/, from December 17, 2004 until February 14, 2014. 

8 The confidence level is defined here as (/ P(x, v)dx, where P(y 2 , v) is the probability density func¬ 


tion of y 2 with V degrees of freedom Press et al. 1992 
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Fig. 2 The duration distribution, ^A/qrb A ^90 of BATSE (blue), Fermi (red) and Swift (Green) GRBs. The 
different curves are shifted in order to avoid overlap. The data bins are evenly spaced in logarithmic scale 
with A log(r 90 ) = 0.1. Bins with less than 5 events are combined with their neighbors in order to achieve 
statistical significance. The black horizontal lines mark the bins that fit a plateau at a confidence interval 
up to 2a. (This figure is taken from |Bromberg et al.|2015) . 


in the BATSE data (20.75/12 x 2 /d.o.f), 0.4-17 s in the Fermi- GBM data (8.7/10 
X 2 /d.o.f) and 0.2-20 s in the Swift data (9.04/8 X 2 /d.o.f). 

Taking a median redshift of z — 2 for Swift GRBs and z. — 1 for Fermi- GBM 
and BATSE bursts, [Bromberg et al.] ( j2015j ) find that in the GRBs’ rest frame, these 
plateaus extend up to intrinsic durations of 7-12 s, consistent with the results ob¬ 


tained by Bromberg et al. 


(2012). Note that the actual t. may be somewhat longer 


than the duration that marks the end of the plateau, but it cannot be shorter. We use 
the duration interval of 7-12 s as our best estimate for the typical t b . 


3.1.3 Implications: the Jet Composition at Early Times 

Eqs. 0, ([8]) and ( | 1 2\ use parameter values inferred from typical GRB jets, after ac¬ 
counting for the jet opening angle (e.g., Bloom et al .|[2003j ). From these equations it 
can be seen that a breakout time of 7-12 s is consistent with that expected for a hy¬ 
drodynamic jet from a WR star with a radius of ~ (3 — 6 )Rq. On the other hand, these 
breakout times are too long for typical parameters expected for a Poynting dominated 
jet. To account for these breakout times, the light cylinder of the compact object at 
the base of the jet should be of the order of r L ~ (2.5 — 5) x 10 8 cm, corresponding to 
an angular frequency of Q nl ~ 60 — 120 rad/s at the base of the jet. Such a frequency 
is too low to allow the engine to power a typical GRB jet ( {Bromberg et al.(2015j l. 
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Fig. 3 The duration distribution, dN(^^/dT )t] of the soft GRBs. The analysis is the same as in Figure 
[ 2 ] only the data from each satellite contains only events that are softer than the median hardness of the 
long GRBs with durations T 90 > 20 s. For BATSE, this corresponds to GRBs having a hardness ratio 
HR32 < 2.6, for Fermi the GRBs have a power law spectral index < —1.5, and for Swift the GRBs have a 
power law spectral index < —1.7. The analysis here updates the analysis in Bromberg et al. 12013} using 
a more complete recent sample. (This figure is taken from Bromberg et al. 2015 1) 


The conclusion arising from this analysis is that during most of its propagation 
within the star the jet has a low magnetization and it propagates as a hydrodynamic 
jet (unless 3D effects significantly increase f bmag ). This result leads to some interest¬ 
ing implications for the properties of long GRB central engines and the conditions 
at the base of the jets. One possibility is that the jet is launched hydrodynamic ally at 
the source. The most probable process for that is neutrino-anti-neutrino annihilation 
above the ro tational axis of the central engine (e.g. Eichler et al.|1989 Levinson and 


Eichler 1993). This scenario can work only if the accretion rate is > 0.1 M 0 s 


so 


that neutrino emission is large enough to power the observed jets (Kawanaka et al. 
|2013[[Levinson ancfG lobus 2013 ). The high accretion rate must be sustained through¬ 
out the entire duration of the GRB, which can last from tens to hundreds of seconds. 


Though a duration of < 30 s seems to be consistent with such a model (e.g. Lindner 
et al. 2010), its seems unlikely to be capable of powering longer duration GRBs. 


A second possibility is that the jet is launched Poynting flux dominated but it dis¬ 
sipates most of its magnetic energy close to the source, and it then propagates as a hy¬ 
drodynamic jet. An appealing process for such efficient dissipation is the kink insta¬ 
bility ( |Lyubarskij|1992| Eichler|1993[ Sprui£etak|1997^ Begelman|1998[ Lyubarskii 


|1999[ |Giannios and Spruit|2006} .' Bromberg et al.| ( 2014} , however, have shown 

analytic considerations, that collapsar jets are less likely to be disrupted by the kink 
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instability. Thus a different process, possibly internal to the jet, may be needed to dis¬ 
sipate the jet energy. A definite answer will be obtained only via 3D numerical sim¬ 
ulations, which are underway. In one such work, Bromberg & Tchekhovskoy (2015, 
in prep.) show that indeed kink instability is unlikely to disrupt a typical collapsar 
jet. Nevertheless, kink modes can grow internally in the jet and lead to efficient dis¬ 
sipation of the magnetic energy via reconnection of the magnetic field lines without 
compromising the jets’s integrity. The outcome of such dissipation is a jet with an 
equipartition between thermal and magnetic energy, which propagates more or less 
like a hydrodynamic jet. 

A third possibility is that the jet changes its character with time. Our conclusion 
concerning the jet composition applies only to the initial phase, while its head is still 
within the stellar envelope. This phase, which lasts ~ 10 s, must be predominantly 
hydrodynamic. Once the jet has breached the star it can be Poynting flux dominated. 
This would require a more complicated central engine that switches from one mode 
to another. While this seems contrived, remarkably, some magnetar models suggest 
such a possibility ( jMetzger et al.|20fT) . One can also imagine accretion disk models 
that initially cool via neutrinos and later on as the accretion rate decreases, become 
Poynting flux dominated ( Kawanaka et al.|2013| >. However, all such models require 
some degree of coincidence as the central engine does not receive any feedback from 
the propagating jet and there is no a priori reason that the transition from one com¬ 
position to the other would take place just at the right stage. 


3.2 Jet steady State Acceleration 


Magnetic acceleration and thermal acceleration are the two main competing mecha¬ 
nisms for the acceleration of GRB jets or outflows. Thus, the acceleration mechanism 
is tightly related to the outflow composition and in particular its degree of magnetiza¬ 
tion, which is both highly uncertain and of great interest. In other sources of relativis¬ 
tic jets or outflows, there are currently better constraints on the composition. Pulsar 
winds are almost certainly Poynting flux dominated near the central source. This 
most likely also holds for active galactic nuclei (AGN) and tidal disruption events, 
as in these sources the central accreting black hole is supermassive, and therefore 
even close to it the Thompson optical depth, Ty, may not be high enough for thermal 
acceleration by radiation pressure (the main competition to magnetic acceleration) to 
work efficiently (e.g., Ghisellini|2012) l. In GRBs or micro-quasars, however, thermal 
acceleration could also work (since Ty 1 is possible, or even likely close enough 
to the source), and the dominant acceleration mechanism is less clear. 

First, let us consider the thermal acceleration of a steady, axisymmetric, and un¬ 
magnetized flow that is initially relativistically hot with p 3> pc 1 . Let the jet cross 
section b eZ«r where r is its cylindrical radius. The relativistic equation of state im¬ 
plies p <x p 4 / 3 , while mass and energy conservation read F pcE = const and F 2 (pc 2 + 
4p)cL = const, respectively (where we have assumed a relativistic velocity, fi = 
v/c ~ 1). The ratio of the two last expressions gives the Bernoulli equation - the 
total energy per unit rest energy (which is conserved without any significant energy 
losses from the system), (1 +4p/pc 2 )F = const. As long as the flow is relativisti- 
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cally hot (p pc 2 ) it accelerates as F p/p p" 1 / 3 oc F 1 / 2 r. This reproduces 
the familiar result for a spherical or conical flow for which F i.e. the Lorentz 

factor grows linearly with the distance z from the central source. Therefore, thermal 
acceleration is relatively fast, efficient and robust. 

Let us now do a similar simple analysis for a cold and initially highly-magnetized 
flow, with (To = Bq/471Pqc 2 1 (e.g., Komissarov 2011). Let the flow be steady. 


axisymmetric, and ideal MHD (i.e. without magnetic dissipation). Let us consider 
the flow between two magnetic flux surfaces defined by r and r + Sr (which are both 
functions of z). Flux freezing (ideal MHD) implies that the poloidal and tangential 
magnetic field components scale as B p 1 /rSr and B p °c 1 /Sr, respectively, in the 
lab frame. Therefore, the tangential field component rapidly dominates far from the 
source, so that B ss Bq sa TB' where B' in the magnetic field in the comoving frame 
of the outflowing plasma (in which the electric field vanishes). Altogether this gives 
B =rB' oc 1 /Sr. Mass and energy conservation read F pcE = const and F 2 (pc 2 + 
B' 2 /4n)cE = const, respectively, where E oc rSr. Their ratio implies a total energy 
per unit rest energy of (1 +cr)F = (1 + (7 o)Fo = fjn ax where cr = B' 2 /Anpc 2 °c r/F Sr 
is the magnetization parameter. Therefore, this results in the following Lorentz factor 
evolution: 


F / (5r 0 r 

fr — 1 + Oq ( 1 
F 0 \ r 0 or 


r 

Fnfix 


= 1 - 1 - 


Sr 0 r 
ro Sr 


(14) 


This immediately implies that for a conical (or spherical) flow, in which Sr °c r and 
Srjr = 5/'()//■(), the Lorentz factor essentially remains constant, F k, const, and the 
flow hardly accelerates. This result can be understood by simple energy considera¬ 
tions. As long as there is no expansion along the direction of motion, the volume of a 
fluid element scales as °c r 2 while its magnetic energy density scales as oc B 2 oc r~ 2 , 
implying a constant magnetic energy and no conversion into kinetic energy. 

More generally, Eq. ( p~4| ) implies that in order for the flow to accelerate, r/Sr must 
decrease, i.e. streamlines must diverge faster than conical. For power-law streamlines, 
z = Zo(r/ro) a = zo[(r + Sr)/{ro + 5ro)]“, one has r/Sr = ro/Sro so there is still no 
acceleration. If one allows the power law index to vary with ro = r(zo), i.e. 0! = afro), 
then one finds Sr/r = (5ro/ro)[l — roa , (ro)o: _2 (ro)ln(z/zo)], and the condition for 
acceleration becomes a' = da/dro < 0. Altogether one can see that such steady-state, 
axisymmetric ideal MHD acceleration is quite delicate and requires a very particular 
configuration of the magnetic field lines. Satisfying this requirement is not trivial, and 
in particular it requires lateral causal contact across the jet. 

A key open question regarding outflows that start out highly magnetized near 
the central source is how they convert most of their initial electromagnetic energy to 
other forms, namely bulk kinetic energy or the energy in the random motions of the 
particles that also produce the radiation we observe from these sources. It is suggested 
by observations of relevant sources, such as AGN, GRBs or pulsar wind nebulae that 
the outflow magnetization is rather low at large distances from the source. This is the 
essence of the well-known a problem - how to transform from cr 1 near the source 
to a <§C 1 very far from the source. 

It has been shown early on that a highly magnetized steady spherical flow acceler- 

1/3 2/3 

ates only up to an asymptotic Lorentz factor F» ~~ (T () , an d magnetization <T.„ ~ cr () 
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( (Goldreich and Julian|l970j ) where Oo 1 is the initial value of the magnetization 
parameter cr, implying that most of the energy remains in electromagnetic form (a 
Poynting flux dominated flow). This is valid for any such unconfined flow, i.e. where 
the external pressure is effectively negligible. A sufficiently large external pressure 


can help collimate and accelerate the flow. It has been found (Lyubarsky 2009 2010b 


Komissarov et al. 


2009 that for a power law external pressure profile, p ext < 


, the 

collimation and acceleration can proceed in two distinct regimes. 

For k>2, the weak confinement regime, the external pressure drops fast enough 
such that the flow loses lateral causal contact while it is still highly magnetized, 
and from that point on it becomes conical and essentially stops accelerating. This 
collimation-induced acceleration increases L x , and decrease ci x , by up to a factor of 
~ 0- ~ compared to the unconfined (quasi-spherical) case, where is the asymp¬ 
totic jet half-opening angle. This arises because lateral causal contact in the jet is 
maintained as long as 0, does not exceed the Mach angle, 0j < 0 M ~ cr’^/r, where 
energy conservation implies aT ~ <7o (for Co 1 and fo ~ 1) as long as the flow 
remains highly magnetized (cr 3> 1). 

For K < 2, the strong confinement regime, the external pressure drops slowly 
enough that the jet maintains lateral causal contact throughout its collimation-induced 
acceleration process. In this case about half of the initial magnetic energy is converted 
into kinetic energy and the flow bacomes only mildly magnetized, cr*, ~ 1, while the 
Lorentz factor approaches its maximal possible value, IT ~ (To. In this regime the 
collimation and acceleration proceed as T r z K ' 4 and the jet remains narrow, 

r»e ) ~ 1 . 

The main problem with this picture, however, is that even under the most fa¬ 
vorable conditions the asymptotic magnetization is (T*, > 1, which does not allow 
efficient energy dissipation in internal shocks within the outflow ( I^ubarekyj2009 


2010b Komissarov et al. 2009| l. It has been found (Tchekh ovskoy et al.|2010 Komis- 
et al.|2010|l that a sudden drop in the external pressure, as may occur when a 


sarov i 


GRB jet exits its progenitor star, can result in a sudden additional acceleration that 
can lead to Fjo0j ^ 1 as inferred in GRBs, but still with o„. > 1. 

These important limitations of the “standard” steady, axisymmetric and non- 
dissipative (or ideal MHD) magnetic acceleration have, on the one hand, led to the 
suggestion that the jets might remain Poynting flux dominated at large distances from 
the source and the observed emission is the result of magnetic reconnection events 
rather than internal shocks (BI and ford|2002[ [Lyuti kov and BI andford|2003] |Lyutikov| 


20061. On the other hand, other models suggested increasing the acceleration effi¬ 


ciency by relaxing one of the standard assumptions, such as axi-symmetry - lead¬ 
ing to non-axisymmetric instabilities that randomize the magnetic field orientation 
(Heinz and Begelman 2000} . Since a highly tangled magnetic field effectively be¬ 
haves like a relativistic fluid (with an adiabatic index of 4/3) this leads to efficient 
acceleration, similar to thermal acceleration of relativistic outflows. What is more, 
both the kink instability mentioned above ( Drenkhahn and Spruit|2002~| , as well as 
other instabilities (such as the Kruskal-Schwarzschild instability in a striped wind; 
Lyubarsky|2010a i can lead to magnetic reconnection, i.e. gradual magnetic dissipa¬ 


tion, which in turn enhances the acceleration due to the conversion of magnetic to 
thermal energy, where the thermal pressure efficiently accelerates the outflow. 
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Fig. 4 Left: Test case for impulsive magnetic acceleration: the energy-weighted mean Lorentz factor (JT) 
of a finite cold shell of plasma initially uniform (with width lg, rest mass density pg and magnetic field 
Bo), highly magnetized (ffo = Bg/4jtpgc 2 2> 1; Co = 30 was used here) and at rest, whose back leans 
against a conducting “wall” while its front faces vacuum, versus the time t in units of the shell’s initial 
fast magnetosonic crossing time tgailg/c. The analytic expectations (dotted and dashed-dotted lines) and 
the results of numerical simulations (diamond symbols joined by a solid line ) are in very good agreement. 
(This Figure is taken from |Granot et al.|201 1 ) . Right: Evolution of the typical (or energy-weighted average) 
Lorentz factor 1 with the distance Raid from the central source, for a finite shell similar to that described 
in the left panel, but for a spherical shell propa gatin g int o an external medium with a power-law density 
profile, p ext = ,4 R k . (This figure is taken from|(jranot|2012a|. 


3.3 Impulsive Magnetic Acceleration 


Replacing the usual steady-state assumption by strong time-dependence is a natu¬ 
ral alternative. This impulsive regime was sparsely studied, and mainly in the non- 
relativistic case (Contopoulos 1995[ >. Recently, a new impulsive magnetic accelera¬ 
tion mechanism was found that operates in the relativistic case (G ranot et al.|201 f) , 
which can be much more efficient than magnetic acceleration in steady flows, and 
can lead to low magnetizations, a <C 1, thus enabling efficient dissipation in internal 
shocks. This qualitatively different behavior of impulsive outflows can be very rele¬ 
vant for GRBs, as well as for other relativistic jet sources such as tidal disruptions or 
flares in AGN or micro-quasars, or even giant flares in soft gamma repeaters (SGRs, 
thought to be magnetars), which also triggered renewed interest in this topic (e.g., 
Levinson|20 1 0| |Lyutikov|20TT Granot|2012a|b[ Komissarov|2012 1 . 

Figure [4] ( left panel) shows results for an impulsive magnetic acceleration test 
case: a cold, initially uniform plasma shell (of with width Iq, rest mass density po 
and magnetic field Bp), highly magnetized (Go = Bo/Atzpoc 2 1) and at rest, with 
a conducting “wall” at its back and vacuum in front of it. A strong, self-similar rar¬ 
efaction wave forms at its front (vacuum interface) and propagates towards its back, 
reaching the wall at t = to « 1q/c. By this time the shell’s energy-weighted mean 


Lorentz factor and magnetization are (F) ~ a,[ ’ and (g) ~ Oq / j . At t > to the shell 
detaches from the wall, keeps an almost constant width (1 ss 21(f) and accelerates as 
(F) ~ Gq/(g) ~ ((Tof/fo) 1 / 3 U P t° the coasting time t c = G^to- At t > t, the shell 
coasts at (F) ~ <7o, its width grows (l/2lo ~ t/t c ) and its magnetization rapidly de¬ 
creases ((a) ~ t c /t), leading to complete conversion of magnetic to kinetic energy 
that allows strong internal shocks to form that can lead to large radiative efficiencies. 


t 2/3 

























GRBs as Sources of Strong Magnetic Fields 


19 


3.4 Interaction with the External Medium and the Reverse Shock 


Let us now consider the evolution of a similar shell in spherical geometry that prop¬ 
agates into an external medium with a power-law density profile, p ex t = ARr k , fol¬ 
lowing Granot ( 2012a] >. The main results are shown in the right panel of Figure [4] 
The initial shell magnetization do and density po °c 1 /do are allowed to vary while 
keeping fixed the values of the initial time or length scale (to ~ Ro/c or Rq), energy 
(E ~ Lto ~ LRq/c or power L), and external density (k <2 in this figure, and A or 
Pext(fio) =ARy k ), which imply fixed r cr ~ (/odo) 1 ^ 8-2 ^ where f 0 = po/p e xt(Ro) 
and R CI ~ /?oLr • Shown are the two dynamical regimes most relevant for GRBs. The 
purple line shows regime 1(1 < do < Lr or a sufficiently low external density) where 
the shell initially expands as if into vacuum (as described in the left panel) and only 
after becoming kinetically dominated and expanding radially is it significantly decel¬ 
erated by the external medium through a strong relativistic reverse shock, which can 
produce a bright emission that peaks on a timescale larger than the duration of the 
prompt GRB emission (the familiar low-d “thin shell”; Sari and Piran|199~5) >. Even¬ 
tually, most of the energy is transferred to the shocked external medium and the flow 
approaches the[Blandford and McKeejlj 1976) self-similar solution. 

3(4— k) 

In regime II (1 < Lr < do < Lx ), depicted by the green line in the right 
panel of Figure [4] the shell is significantly affected by the external medium while 
it is still Poynting dominated (at R > R u ~ Ro(fo<J 0 I ^ 3 ) 3 /( 10—3 *)), thus suppress¬ 
ing the reverse shock (which is either non-existent or very weak). The shell remains 
highly magnetized and gradually transfers its energy to the shocked external medium 
through pdV work across the contact discontinuity up to R cr , after which the flow ap¬ 
proaches the Blandford-McKee solution. In this regime no significant reverse shock 
emission is expected, and the onset of the afterglow (i.e. the peak of the emission 
from the shocked external medium) is expected to be on a timescale comparable to 
the prompt GRB duration (i.e. a high-d “thick shell”). 


In addition, there are other regimes not shown in this figure, but all of the regimes 
are mapped in the relevant parameter space in Figure]^] In regime III (1 < Lx' 4 k ^~ < 
do) the external density is high enough that there is no impulsive acceleration stage 
where (T) R 1 / 3 , and instead (T) ~ do/(d) R( k ~ 2 )/ 4 at Ro < R < R c r ~ Rdec> 

and then approaches the Blandford-McKee solution (its observational signatures are 
expected to be similar to regime II). In regime IV (Lr < 1) the external density is so 
high that the flow remains Newtonian all along (as might happen while the GRB jet 
is propagating inside a massive star progenitor). There is also an “exotic” regime II* 
that exists only in a highly stratified external medium (10/3 < k < 4). 


Under realistic conditions, GRB variability times are in practice typically large 
enough that the flow should first undergo quasi-steady collimation-induced acceler¬ 
ation that saturates, and only later the impulsive acceleration kicks in and operates 
until the flow becomes kinetically dominated (see Figure |6j. Moreover, one typically 
expects the outflow from the central source to consist of many sub-shells rather than 
a single continuous shell. The effects of such multiple sub-shells in the outflow can 
be important, and the collisions between them may provide efficient energy dissipa¬ 
tion that can power the GRB emission (|Granot|2012b[ Komissarov|2012|l. They may 
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Fig. 5 Phase space diagrams of the different dynamical regimes: in the fa (To plane for k < 10/3 (top left 
panel), F c r -Gb plane for k < 10/3 (top right panel), fo~r CT plane for k < 10/3 (bottom left panel), and in 
the fa do plane for 10/3 <k <4 (bottom right panel). Each regime is labeled and denoted by a different 
color, and the borders between the different regimes are indicated by labeled thick black lines. (This figure 
is taken from |Granot|20 1 2a| ) 


also allow a low-cr “thick shell”, i.e. a strong relativistic reverse shock peaking on a 
timescale comparable to the prompt GRB emission, which is not possible for a single 
shell. For a long-lived source (e.g. AGN) with initial sub-shell widths U) and separa¬ 
tions /g a p, each sub-shell can expand by a factor of 1 + Zg a p//o- Its magnetic energy 
decreases by the same factor (where a «, ~ /o/ZgapX and may be converted to kinetic 
or internal energy, or radiation. For a finite source activity time, the merged shell 
can still expand further and convert more magnetic energy into other forms (even 
without interaction with an external medium). Important related points that warrant 
further study are the transition from quasi-steady collimation-induced acceleration to 
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Fig. 6 An illustration of the expected transition from (quasi-) steady collimation-induced acceleration 
near the central source (in red ) to impulsive acceleration further away from the source (in green). The 
energy weighted mean Lorentz factor F is shown against the distance R from the central source, and a few 
critical radii and Lorentz factors are indicated following the notations of Granot et al. (2011) and |Granot| 
(20l2a) . 


impulsive acceleration, both in a single shell and in multiple sub-shells, as well as the 
dissipation in the interaction between sub-shells and its effect on the outflow acceler¬ 
ation and the resulting emission, such as a possible photospheric spectral component. 


4 Dissipation and Prompt Emission 


As discussed above, GRBs must be associated with relativistic outflows ejected by 
a stellar mass compact source, with a bulk Lorentz factor F > 100 in order to avoid 
the compactness problem (Baring and Harding 1997 1 Lithwick and Sari 2001; Granot 


et al|2008| jHa scoet et al.|2012) . This also naturally explains the afterglow through 

the deceleration of the ejecta by the external medium, whereas the observed fast 
prompt variability implies that the prompt emission must be produced by internal 
dissipation within the ejecta (Sari and Piran 1997) ). Therefore, the analysis of the 
GRB prompt emission may provide valuable information on magnetic fields within 
an ultra-relativistic jet. It can put unique constraints on the state of the jet at the end 
of the acceleration phase, and more specifically on the geometry of the magnetic field 
and the magnetization at a large distance to the central source, where the y-ray emis¬ 
sion is produced. This is, however, a difficult task as it requires a full understanding 
of the nature of the dissipative mechanisms and of the radiative processes at work. 

There are several possible emission sites for the GRB prompt emission: (i) a com¬ 
ponent can be emitted at the photosphere, where the ejecta becomes transparent to its 
own photons; another component can be produced above the photosphere in the op¬ 
tically thin regime, associated with either (ii) internal shocks propagating within the 
ejecta (Rees and Meszaros 1994); or (iii) magnetic reconnection (Thompson 1994; 
Spruit eTal.|2001| l. These three dissipative mechanisms extract energy, respectively, 
from the thermal, kinetic or magnetic reservoirs. The expected prompt emission com¬ 
ponents are therefore strongly related to the composition of the ejecta. 
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4.1 Photospheric Emission 

4.1.1 Noil-Dissipative Photospheres 

When internal dissipation below the photosphere is negligible, the expected thermal 
emission at the photosphere is well understood with precise predictions ( jPaczynski 
1986 Goodman|1986| Shemi and Piran|1990 Meszaros et al.|1993) . Only few the 


oretical uncertainties remain, mainly related to the lateral structure of the jet (see 
e.g. Lundman et al.|20l3| Deng and Zhang|2014| l. Assuming that the photosphere is 
above the saturation radius, the photospheric radius is given by (Meszaros and Rees 


2000 Meszaros et al.|2002] Daigne and Mochkovitch|2002||Hascoet et al.|201 3]i 

KE _ .„i3 Kfl.2^53 


^ph — 


87TC 3 r 3 (l 


~3x 10 
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cm 


(15) 


where K 0.2 is the matter opacity in units of 0.2cm 2 g 1 , T = IOOI 2 is its bulk Lorentz 
factor, E = 10 53 Es3 erg s 1 is the isotropic equivalent jet power, and a is the mag¬ 
netization parameter at the end of the acceleration, so that E/{\ +a) is the isotropic 
equivalent kinetic power. The observed photospheric luminosity and temperature are 

\" 2 / 3 T „ / tf . X - 2 / 3 


f'ph — ^th E 


R 


ph 


rRo) 


and 


Tph — 


To (R p h\ 
1+z V ERo) 
1/4 


(16) 


where the initial temperature equals To = (ethE/A tioRq) . Here e t h is the thermal 
fraction of the jet power at the base of the flow, located at the initial radius Rq, and 
1 — £,h is therefore the initial magnetic fraction. In the case of a passive magnetic field 
carried by the outflow without contributing to the acceleration, the initial thermal 
fraction £ t h and the magnetization at the photosphere cr are related by (7 pass i ve = Oq = 
(1 — £th)/£th- An efficient magnetic acceleration leads to cr < CT pass j ve (Spruit et al. 


|2001| ). The predicted spectrum is quasi-thermal, with an exponential cutoff at high 
energy and a power law at low-energy with a photon-index of a ~ 0.4, which differs 
from the a = 1 slope of the Raleigh-Jeans spectrum due to the peculiar geometric 
shape of a relativistic photosphere (Goodman 1986; Beloborodov|201 1), 


4.1.2 Dissipative Photospheres 

If dissipation occurs below the photosphere, the emitted spectrum can be significantly 
different than the previous case: a high-energy tail can be produced by comptoniza- 


tion due to the presence of relativistic electrons ( Thompson 

1994; 

Meszaros and Rees 

2000 Rees and Meszaros 2005 Giannios and Spruit 2007 

Belo 

Gorodov 2010), and 

the low-energy slope can be modified by synchrotron radiation (Pe’er et ak|2006 


|Vurm et ~aI7||2011| >. The resulting observed spectrum may now appear non-thermal, 
with several components. This scenario is more uncertain than the previous one. The 
nature of the sub-photospheric dissipative mechanism must be identified, with several 
candidates: early internal shocks, gradual magnetic reconnection (Thompson! 1994 


Giannios and Spruit 2007), neutron-proton collisions (Beloborodov||2010|l, etc. An 


important relevant issue is related to the photon production efficiency and thermal- 
ization deep within the ejecta (Vurm et al .j20i3j l. 
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4.2 Non-Thermal Emission in the Optically Thin Regime 
4.2.1 Electron Acceleration and Synchrotron Radiation 


Non-thermal emission can be produced above the photosphere if some internal dissi¬ 
pation processes can lead to efficient electron acceleration. In this case, two natural 
candidates for the dominant radiative process are the synchrotron radiation and the 
inverse Compton scatterings of synchrotron photons by relativistic electrons (SSC). 
However, the measurement of the prompt y-ray spectrum over a broad spectral range 
(keV-GeV) in a few bursts by Fenw/GBM+LAT can rule out the possibility of SSC 
being dominant in the soft y-ray range, as it would lead either to a strong synchrotron 
peak at lower energy, or a strong second inverse Compton peak at higher energy, 
which are not observed (Bos njak et al.|2009 Pi ran et aT| 2009 ). Therefore, the dis¬ 
cussion is focussed on the synchrotron radiation of relativistic electrons, with several 
possible dissipation mechanisms responsible for the acceleration of electrons. 


4.2.2 Internal Shocks 


If the magnetization at a large distance to the central source is sufficiently low, strong 
internal shocks are expected to form and propagate within the ejecta due to the vari¬ 
ability of the ejected outflow. A large range of radii is expected. 


R, 


- 2r 2 c{t mI ->■ f C RB) - ri ( 6 X 1()12 Gar, -2 6 X 10 15 f G RB.l) ™, (17) 


where Aar ,-2 is the shortest timescale of variability, in units of 10 2 s and /grim is 
the total duration of the relativistic ejection, in units of 10s. The shocks are expected 
to be mildly relativistic, except for a very large amplitude of variation of the initial 
Lorentz factor. The dynamics of the internal shocks phase has been studied in detail 
( Kobayashi et al.|1997| Daigne and Mochkovitch 1998 2000), from a simple ballistic 
approximation to a full hydrodynamical code, and is well understood. Up to fa ~ 40% 
of the kinetic energy can be dissipated (for a low magnetization outflow, cr -C 1), 
depending again on the initial distribution of the Lorentz factor. 

On the other hand, large uncertainties on the emission remain, due to the poor 
understanding of the microphysics of mildly relativistic shocks (for a recent review 
of relativistic collisionless shocks see Sironi et al. [2015| . It is usually parametrized 
by assuming that a fraction fig of the internal energy is injected into an amplified 
random magnetic field at the shock, whose structure is not known, and a fraction £ e 
is injected into a fraction £ of electrons, which are therefore accelerated into a non- 


thermal distribution with slope p ( dN e /dy e °< y e l '). To reach the soft y-ray domain 
by synchrotron radiation, the fraction of accelerated electrons must be low, ^ < 10 


,-2 


(Daigne and Mochkovitch|1998| Bosnjak et al.|2009[ Daigne et al.|201 l]|Beniamini 


and Piran 2013). On the other hand, values of £ e close to equipartition (£ e = 1 /3) are 


required to explain the huge luminosities of GRBs. If the radiative efficiency is high, 
a fraction /is ~ fae e ~ 0.01 — 0.1 of the initial kinetic power can be converted into 
radiation in internal shocks. With such assumptions, the non-thermal emission in the 
comoving frame of the shocked regions can be computed with a detailed radiative 
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model including all the relevant processes, namely synchrotron radiation and slef- 
absorption, inverse Compton scatterings and photon-photon annihilation. The contri¬ 
butions of each internal shock can then be added with an integration over equal-arrival 
time surface of photons to the observer in the source frame to produce synthetic light 
curves and spectra that can be directly compared to observations. The predicted spec¬ 
trum shows several components, a strong synchrotron peak in the soft y-ray range 
and a weaker inverse Compton peak at higher energy (|Bosnjak et al.|2009)>. 


4.2.3 Magnetic Reconnection 


If the magnetization at large distances remains high (a > 1), then internal shocks 
are either significantly suppressed or in some cases cannot form altogether (Mimica 


and Aloy 2010; Narayan et al. 2011). In such cases, electrons may be accelerated 
predominantly in magnetic reconnection sites. This scenario is even more uncertain 
(less understood) than the two previous ones (photosphere, internal shocks), but it is 
under extensive investigation (for a recent review on relativistic magnetic reconnec- 
Kagan et al.|2015). Some authors considere d a gradual reconnection starting 


tion see 


at a small radius and extending up to R rec ~ 10 13 cm (Drenkhahn and Spruit 


2002 


Gi- 


annios|2008| ). In this case, most of the dissipation occurs below the photosphere, cor¬ 
responding to the dissipative photosphere scenario discussed above. If, on the other 
hand, reconnection remains inefficient below the photosphere, it could occur at larger 
radii i n the optically thin regime. In the simulations by McKinney and Uzdensky 


(20121), a catastrophic dissipation of the magnetic field occurs at R l ec ~ 10 13 — 10 14 cm 




when reconnection enters a rapid collisionless mode. Zhang and Yan ( 2011J ) pro¬ 
posed another scenario where reconnection is triggered by internal shocks, the so- 
called ICMART model. The typical radius may be as large as R rec ~ 10 15 cm. The 
microphysics in the reconnection sites is also uncertain. One expects many electron 
acceleration sites, which may move relativistically in the outflow’s rest frame. The 
non-thermal electron distribution may be somewhat harder than in shock accelera¬ 
tion (see e.g. Sironi and Spitkovsky: 2014] Kagan et al. 2015] , and the acceleration 
process slower, which can lead to a different shape for the synchrotron emission. Con¬ 
trary to internal shocks, detailed calculations of the light curves and spectra based on 
a detailed radiative model coupled to a dynamical simulation are not yet available. 


4.3 Magnetic Field in Emission Sites 


The models discussed above have very different implications for the magnetic field: 

- Case 1; most of the prompt emission is due to a dissipative photosphere. Then, the 
magnetic field must be large enough at the photosphere to produce synchrotron 
radiation and affect the low-energy spectrum. If this magnetic field is generated 
by the dissipation process (e.g. shocks; Sironi et al. |2015 1, it is most probably 
random. Otherwise, an ordered field must be present. The corresponding initial 
magnetization must either be low (otherwise the photospheric emission is weak), 
or high with very efficient reconnection below the photosphere, which then leads 
to a possible candidate for the sub-photospheric dissipation process. 
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- Case 2: the prompt emission is mostly non-thermal, from an optically thin region: 

- For internal shocks to be the dominant dissipation process, the magnetization 
at large distances from the source must be low. A random field is generated 
locally at the shock front, where the electrons are accelerated. However, the 
magnetic field felt by the radiating electrons must be considered far behind 
the shock front (as the radiative cooling length is much larger than the plasma 
skin depth) where its strength and structure are not well known. 

- If reconnection dominates the dissipation then cr must be large far from the 
source. The ordered field is destroyed at the reconnection sites, but if electrons 
have enough time to migrate from their acceleration site before radiating, their 
emission may still be mostly in the large-scale ordered field. 

Observations of the GRB prompt emission, discussed in the next subsection, can put 
strong constraints on these various emission models. 


4.4 Constraints from the Observed Prompt Soft y-ray Emission 
4.4.1 Light Curves 


All the scenarios discussed above can reproduce the observed variable light curves. 

There are, however, important differences: 

- (Dissipative) photosphere: the emission radius is low (see Eq. Therefore the 
curvature effect, i.e. the spreading of photon arrival times from different angles 
with respect to the line of sight over an angular timescale tg = R/2r 2 c of a flash 
of photons emitted at the same time and radius, is negligible: the observed light 
curve directly traces the activity of the central engine. 

- Internal shocks: the light curves trace the source activity (Kobayashi et al. [ 1997) 
Daigne and Mochkovitch |1998), but two effects now affect the observed pulse 


shapes: the curvature effect (due to a larger radius) dominates the pulse decay 
(Genet and Granot 2009; Willingal e et al.|2010| >, and the radial or hydrodynamic 
timescale due to shock propagation, t r = AR/2r 2 c, dominates the pulse rise and 
overall shape (Daig ne~and Mochkovitch| 1 998]|2003[ Bosnjak and Daigne|2014| l. 

- Reconnection: again, the light curve traces the source activity, with new effects 
due to relativistic bulk motion in the local jet’s frame. Relativistic motions of 
emitting plasma in the jet’s frame cause rapid variability (that should show up as a 
distinct component in the Fourier power spectrum), while a slower envelope may 
arise from their combined effect (Zhang and Zhang |20l4] ) or from slower emitting 
plasma. This can be tested by characterizing the observed variability. Analysis of 


GRB light curves shows a continuum of timescales (see e.g. Beloborodov et al. 

2012[>, which does not support the reconnection model of 


2000 Guidorzi et al. 


Zhang and Zhang (2014) (see however Gao et al.|2012[). A possible concern ap¬ 


pears if the emission is produced by many relativistically moving emitters: the 
predicted pulse shape may be too symmetric compared to observations (jLazar 
et al.[2009]). However, both concerns (the power spectrum and pulse shapes) may 


be solved if the reconnection occurs in relatively ordered thin layers located be¬ 
tween anti-parallel regions in the outflow (with a geometry of thin quasi-spherical 
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shells) and the relativistic motions in the jet’s frame are limited to these layers 
(Beniamini & Granot 2015, in prep.). Such a model may also account for many 
of the correlations that are observed in the prompt emission. 


4.4.2 Polarization 


Measuring the polarization in the y-ray domain remains challenging. A very large 
degree of polarization was claimed by Coburn and Boggs ( 2003| , but it was later 
refuted by others (Rutledge and Fox |2004 Wigger et al.jl 2004j > as not being statis¬ 
tically significant. Only a few later measurements (by INTEGRAL, GAP) are avail¬ 
able, however with a low or moderate significance (McGlynn et al.|2007 Gotz et al. 


2009; Yon etoku et al.|20lT||2012[|Gotz et al.|2014| l. Such measurements (if reliable) 

can put constraints on the magnetic field geometry in the emission sites (e.g., Gra- 


not and Konigl|2003 Granot|2003[ Lyutikov et al.|2003 Nakar et al.|2003| l. Current 
observations seem to favor synchrotron radiation in an ordered field with patches, 
which would favor emission in the optically thin regime above the photosphere. It is 
however not trivial to justify a highly ordered field in the internal shocks model (a tur¬ 
bulent field at the shock is required for particle acceleration; The structure of the field 
on intermediate scales between the plasma and the dynamical scales is less known), 
or in the reconnection model (in principle the structured field is destroyed by recon¬ 
nection, but the remaining field can still possess significant structure and electrons 
may also radiate somewhat outside of these localized reconnection regions). Either 
way, more definitive polarization observations are needed before strong conclusions 
can be drawn. 


4.4.3 Spectrum 


The prompt soft y-ray spectrum is usually fitted by a phenomenological model in¬ 
troduced by Band et al. ( 1993) >, which consists of two power laws with low- and 
high-energy photon indices of a and /j, smoothly connected at the peak energy 
This eliminates non-dissipative photospheres for the dominant contribution to the 


emission, which puts a strong constraint on the initial magnetization (low Daigne 


and Mochkovitch|2002|l, favoring magnetic acceleration of the outflow. The remain¬ 


ing scenarios for the prompt GRB emission are either a dissipative photosphere, or 
a combination of a weak photospheric emission and a non-thermal component due 
to shocks or reconnection. The discussion is then focussed on the general shape of 
the spectrum, and the low-energy photon index a, which is observed to be close to 
a ~ — 1 ( |Preece et al.|2000||Kaneko et al,|2006l|Nava et al.|201 l~||Gruber et al.|2014] i. 


- Dissipative photospheres: the value of a can be reproduced by adjusting the 
magnetization, which controls the synchrotron emission at low energies. The the¬ 
oretical instantaneous spectral peak is narrower than the observed time-integrated 
spectral peak, but the comparison should be made using a theoretical time-integra¬ 
ted spectrum, which should broaden it. 

- Internal shocks; to reproduce the high luminosities and the short timescale vari¬ 
ability of GRBs, the radiating electrons must be in the fast cooling regime (Cohen 
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et al.|1997] |Sari et al.|1998[ |Ghisellini et al.|2000| l, i.e. their radiative timescale 
must be shorter than the dynamical timescale that governs the adiabatic cooling 
for the spherical expansion. This leads to a predicted photon index a < —3/2, in 
contradiction with observations (the so-called “synchrotron linc-of-deatlf', |Preece 
et al. 1998} . Another potential problem is that the resulting spectrum is too broad 
around the peak. The two problems are naturally connected. Several possibilities 
have been discussed to solve this issue: (i) inverse Compton scatterings in the 
Klein-Nishina regime affect the cooling of electrons, leading to photon indices 


a < — 1 ( Derishev et al.|200T| Wang et al.|2009| Nakar et al.|2009 Daigne et al. 
201 1} . This puts a constrai nt on the strength of die magnetic field, which sh ould 
be small, with £b < 10 3 ( |Daigne et ah 2011 Barniol Duran et al.||2012 i; (ii) 
in the marginally fast cooling regime ( Daigne et al. ||20 1 1| |Beniamini and Piran 
2013} , where the radiative timescale is close to the dynamical timescale but still 
below, the electron radiative efficiency can remain large enough (> 50%) to ex¬ 
plain the observed luminosities, but the synchrotron spectrum is strongly affected: 
the intermediate region of the spectrum below the peak with a photon index —3/2 
disappears and the slope a = —2/3 usually associated with the inefficient slow 
cooling regime is measured. This regime is also favored by weak magnetic fields; 
(iii) in the fast cooling regime, electrons radiate on timescales which are long 
compared to the plasma scale at the shock front, but small compared to the dy¬ 
namical timescale. Then they experience a magnetic held that is not necessarily 
the same as the turbulent held just behind the shock (see the structure of the mag¬ 
netic held in the simulations, e.g. Keshet et al. |2009} . If the held is decaying on 
this intermediate scale, it will affect the synchrotron spectrum and can lead to a 
hard spectrum, with — 1 < CC < —2/3 (|Derishev|2007[ Lemoine|2013||Uhm and 
Zhang|20l4l|Zhao et al.|2014} . 

- Reconnection: the slow electron heating in the turbulent held can lead to hard 
synchrotron spectra with a ~ — 1 ( |Uhm and Zhang||2014} . It is unclear if the 
expected hard power-law index p < 1.5 of the non-thermal electron distribution 
(e.g., Sironi and Spitkovsky |2014| ) can be identihed in the observed spectrum. A 
potential issue is that the presence of many emitting regions that move relativis- 
tically in random directions in the jet’s fame might lead to much broader spectra 
than observed. This may be alleviated, however, in models where these regions 
move predominantly along the thin reconnection layer that is located between re¬ 
gions of oppositely-directed magnetic held in the How, normal to the jet’s bulk 
motion (Beniamini & Granot 2015, in prep.). 

Recently the description of observed GRB spectra in the soft y-ray range has been 
greatly improved by Fermi /GBM observations. An important result is the identifica¬ 
tion of significant deviations from the Band spectrum, which seem to be related to 
the presence of a weak thermal component below the dominant non-thermal one 


(see e.g. Ryde^e^c^[20K) 2011 Gumecjuld j201 f Burgess^etak |20 1 1} | Axelsson| 


et al.|201 2| Guiriec et al.|2013| Burgess et al.|2014 Guiriec et al.|2015} , as illustrated 

in Figure [7] A natural explanation is to associate the quasi-thermal weak component 
to a (non-dissipative) photosphere and the Band component to synchrotron radiation 


9 


Or possibly even a dominant photospheric component in the case of GRB 090902B. 
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Fig. 7 Two examples of quasi-thermal components detected in GRB prompt spectra. Left: a weak quasi- 
thermal component at T ~ 38 keV in the long GRB 100724B ( fro m|Guiriec et al.|2011) . Right: a quasi- 
thermal component at T ~ 12 keV in the short GRB 120323A (from |Guiriec et al.|2013) . 






Fig. 8 An example of a synthetic GRB with the contribution from the photosphere and internal shocks 
computed self-consistently. The photospheric emission is plotted in red, the non-thermal emission from 
internal shocks in blue, the total in black. Top-left: initial distribution of the outflow Lorentz factor at the 
end of the acceleration phase. Top-right: light-curves in the GBM energy range. Bottom-left: spectrum. 
Bottom-right: spectral evolution (observed peak energy of the non-thermal component and temperature of 
the photosphere). The dashed and dotted lines show the expected result when integrating the spectrum over 
different timescales. The parameters are E — 10 53 erg/s, £ t h — 0.03 (high initial magnetization, Gq — 32.3), 
G = 0.1 (low magnetization at large distance), Rq — 3 x 10 7 cm. (All panels are taken from |Hascoet et akj 
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from electrons accelerated either in shocks or in reconnection: Figure[8]shows an ex¬ 
ample of a synthetic burst with these two contributions in the case of internal shocks. 
The weakness of the photospheric emission puts interesting constraints on the initial 
magnetization of the outflow (Daigne and Mochkovitch 2002 Hascoet et a l.|2013| >, 
favoring an efficient magnetic acceleration, with a large range of initial magnetiza¬ 
tion in the GRB population, £ t h < 0.01 (Go > 100) in most cases where no detection is 
made and £ t h — 0.01 —0.1 (ofo — 10— 100) in less frequent cases like GRB 100724B 
(Hascoet et al. 2013). GRB 090902B with £* ~ 0.3 — 1 (<Jq < 2.3) remains an excep¬ 
tion within long GRBs, and the short GRB 120323A appears as an intermediate case 
between GRBs 100724B and 090902B with £ th ~ 0.1 - 0.5 (ct 0 ~ 1 - 9) dGuiriec] 
et al.[2013). The fact that the photospheric emission seems brighter in the only case 


of detection in a short GRB (GRB 120323A, |Guiriec et al ,|2013] ) may indicate a dif¬ 
ferent acceleration mechanism. If this interpretation is correct, these recent detections 
rule out purely thermal acceleration (standard fireball) at least in long GRBs. 


4.5 Constraints on a Poynting Flux Dominated Outflow 


Several authors (Lyutikov 2006. Giannios and Spruit 2006} |Zhang and Yan 20111 
have proposed that GRB jets are Poynting flux dominated all the way up to the 
emission region. The prompt y-rays arise, in this case, from a process that converts 
this magnetic energy to radiation. Obviously, this cannot take place directly and one 
has to invoke some sort of magnetic dissipation (e.g. reconnection) that converts the 
magnetic energy to accelerated electrons (or electron-positron pairs) that emit the 
observed y-rays. Particular support for this idea came with the claim of of strong 
polarization in the prompt emission by Coburn and Boggs (2003), which was later 
refuted (Rutled ge and Fox|2004| |Wigger et al.|2004| l. Such polarization could arise 
if the magnetic field is ordered and this will arise naturally if the magnetic field is 
dominant ( jGranot and Konigl|2003||Granot|2003}|Lyutikov et al.|2003| l. 

However, the efficiency of the synchrotron emission process poses serious con¬ 
straints on models in which the emission region is Poynting flux dominated (Beni- 
amini and Pira n|2014} . Consider a Poynting flux dominated outflow and an observed 
(isotropic equivalent) y-ray luminosity Ly. This luminosity immediately sets a lower 
limit on the strength of the magnetic field B in the rest framep^of the central source, 
Ly < R 2 m B 2 c , where R cm is the emission radius. Accelerated electrons effectively 
emit synchrotron radiation. The critical issue here is that synchrotron emission is too 
efficient. The accelerated electrons cool so rapidly in a strong magnetic field that their 
lower bands (X-rays and optical) synchrotron emission would produce a signal that 
is much stronger than the observed emission in these bands. 

The observed prompt upper limits in the optical or the X-rays set strong con¬ 
straints on the conditions within the emitting region. First, if the observed y-rays are 
due to some other (non-synchrotron) emission process then this process must be ex¬ 
tremely efficient and its cooling time should be significantly shorter than the relevant 
synchrotron cooling time (see e.g. Figure [9]). 


10 The magnetic field in the jet’s frame is B' = B/F, where / is its bulk Lorentz factor. 
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Fig. 9 If the dominant y-ray emission mechanism is not synchrotron, then in order for it to be able to 
tap a significant fraction of the electrons’ energy its cooling time, t c , must be shorter than that due to 
synchrotron radiation. t c ,syn. which is depicted here by the contour lines. The observational constraints 
FV,syn,opt < 1 mJy, F v , syni x- ray < 1 mjy and V syn ,LATF v ,syn,LAT < 10- 7 erg s _1 cm -2 further constrain t c . 
Beyond the corresponding lines, t c should be significantly shorter than f c syn in order for the synchrotron 
not to overproduce the upper limits on the optical. X-ray or GeV fluxes. Within the black region the 
synchrotron emission produces the observed prompt y-rays. The conditions Ty < 1, 10 (dot-dashed lines; 
a Thompson optical depth that is not too large) define general limits on the parameter space. (see |Beniamini| 
|and Piran|2014| for more details). 


Alternatively, if the observed prompt y-ray emission is synchrotron then there 
must be a rapid reaccelerating process that keeps the electrons with the right Lorentz 
factor so that they would not cool too much and emit strongly in lower energy bands, 
in particular in soft X-rays ( |Ghisellini and Celotti|1999[|Kumar and McMahon|2008| 
Fan |201(i| i. This requires strong fine tuning as the Lorentz factor range in which the 
electrons must be kept is rather narrow (a factor of ~ 3-10) . Multi-zone configu¬ 
rations in which the electrons escape the emitting region before cooling and over¬ 
producing X-ray or optical emission are also a possibility (e.g. as in the ICMRAT 
model [Zhang and^Ya n 2011). Beniamini an d Piran| ( |2014| ) considered several such 
two-zone toy models (in which electrons are accelerated in one region and emit in 
the other) but proper conditions could not be found in any of them. 

These considerations pose severe constraints on prompt emission models that in¬ 
volve Poynting flux dominated outflows. Any emission model in such a regime should 
satisfy these constraints. Lacking a model that satisfies all these constraints, it is likely 
that if the outflow is initially Poynting flux dominated then the magnetic energy is 
dissipated before the emitting region, where it must be subdominant. 
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4.6 Spectral Diversity - Spectral Evolution 


The peak energy E v varies a lot from one GRB to another, from a few to tens of keV 
(X-ray Flashes, X-ray Rich GRBs; [Sakamoto et al.|2005j ) to over 10 MeV ( |Axelsson| 
et al.||2012). An important property is that short GRBs are harder with larger peak 


energies (Kouveliotou et al. 1993 Guiriec et al. |2010j i. Spectral evolution is also 
found: when time-resolved spectroscopy is possible, the GRB spectrum is always 
found to strongly evolve during the prompt phase (see e.g. Lu et al. 2012 Burgess 


et al.|2014[ Guiriec et al. [201 5] for recent analyses of Fermi GRBs, or Preece et al. 


2014 for a very bright case where the spectral evolution can be studied in great detail); 


£p typically varies over more than a factor of 30 within an individual pulse. Spectral 
and temporal properties appear correlated within GRB pulses: hardness following the 
intensity, pulses being narrower and peaking earlier at higher energies, etc. 

Both the spectral diversity between different GRBs and the spectral evolution 
within individual GRBs are hard to reproduce by models. Reconnection models are 
barely developed enough to allow discussion of these observations. In dissipative 
photospheres, variations in E v are related to changes in the properties of the out¬ 
flow ejection leading to a change in the location of the photosphere (|Pe’er|2008 Be¬ 


loborodov 2013; Deng and Zhang 2014). A potential issue is to explain how the dis¬ 
sipative process adjusts to always remain located just below the photosphere (unless 
it always occurs over a wide range of radii, in which case it should also occur above 
the photosphere, so this would not be a pure photospheric model). In internal shocks. 


the spectral evolution is reproduced qualitatively (Daigne and Mochkovitch 1998 
2003[|Asano and Meszaros |201 1| Bosnjak and Daigne| 2014 1 , and can even be repro¬ 

duced quantitatively with some constraints on microphysics parameters (Bosnjak and 


Daigne 2014), which may indicate non-universal values in mildly relativistic shocks, 


as : 


suggested for instance by Bykov and Meszaros ( 1996| see also Bykov et al.|2012|l. 


The spectral diversity is also naturally explained by variations in the lifetime and 


variability of the central engine ^Barraud et al. 2005 

). The hardness-duration relation 

is well reproduced (Daigne and Mochkovitch 1998 

Bosnjak and Daigne 2014). 


4.7 Constraints at Other Wavelengths 

The discussion above was centered on observations in the soft y-ray range, where the 
prompt emission is observed in most GRBs. We discuss here briefly some additional 
constraints coming from observations of the prompt emission at other wavelengths. 


4.7.1 The End of the Prompt Emission: the X-ray Early Steep Decay 


Swift/X RT discovered in most GRB X-ray afterglows an early steep decay at the end 
of the prompt phase, before recovering a plateau and/or a standard afterglow decay 
( Nousek et al.|2006[|Zhang et al. 2006 O’Brien et al.|2 006 ). A natural explanation is 
provided by the high-latitude tail of the prompt emission, once the on-axis emission 
has stopped (Kumar and Panaitescu|2000|[Genet and Granot|2009jl. It can reproduce 
the observed temporal decay and spectral evolution (Liang et ak[|2006[ Willingale 
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et al.] [2010[ ). This puts a strong constraint on the emission radius at the end of the 
prompt phase, which may be fulfilled by internal shocks, and also possibly by recon¬ 
nection models if the radius is large enough (Hascoet et al. 2012) . On the other hand, 
it is incompatible with photospheric models, which must instead explain the rapid 
decay phase by a universal behavior of the central engine when it is switching off. 


4.7.2 Prompt GeVEmission 


FermiTLKV detects GeV emission in some GRBs ( |Ackermann et al.||2013| >. As de¬ 
tection requires enough photons in its energy range (tens of MeV to > 300 GeV), it 
detects mainly very bright GRBs, in terms of both their GeV fluence and their total 
fluence (and thus also in terms of Fy ls() ). For the same reason, LAT detects a smaller 
fraction of short GRBs compared to soft y-ray instruments, since their fluence is 
typically much smaller than that of long GRBs. Bright enough LAT GRBs show a 


distinct high-energy spectral component, usually fitted by a power law (e.g., Abdo 
et al.||2009t Ackermann et al.||2010 2013| >. The observed variability in the prompt 
LAT light curve indicates an internal origin. It is followed by a long-lasting emission 
(with a power law in time and energy) that likely originates from the deceleration 
phase or early afterglow. In dissipative photospheric models, it is hard to produce 
GeV photons due to strong yy annihilation. However, additional processes such as 
later scatterings of prompt photons by the external medium can explain this GeV 
emission (see e.g. |Beloborodov et al.|2014] ). In reconnection models, spectral models 
cannot make such predictions yet. In internal shocks, such multi-component spectra 
are expected (Guetta and Granot 2003; Bosnjak et al. 2009; Asano and Meszaros 
201 2[ Bosnjak and Daig ne|2014 l; the fact that the GeV component is usually weaker 
than the soft y-ray component co nstrains the strength of the magnetic field, imp lying 
that it must be weak, £b < 10 2 (Daigne et al. 2011 Bosnjak and Daigne 2014 1 . 


4.8 Prompt Emission Summary 

The dissipation mechanism and radiative processes responsible for the prompt GRB 
emission are still not well understood due to the complex physics involved, both 
on large and micro-scales. The lack of strong thermal components in GRB spectra 
suggests a high initial magnetization in GRB outflows, while prompt GRB observa¬ 
tional constraints imply a low magnetization in the emission region. Put together, this 
strongly suggests either very efficient conversion of magnetic to kinetic energy, which 
leaves a low magnetization in the emission region, and allows for efficient internal 
shocks ( Granot etal.|20 11 Granot 201 2bj i, or strong magnetic reconnection that con¬ 
verts magnetic energy to thermal energy and accelerates particles, and yet somehow 
leaves a low enough magnetization where these particle radiate most of their energy. 

A weak thermal emission can be produced at the photosphere, followed by a 
dominant non-thermal emission at larger radii. Depending on the efficiency of the ac¬ 
celeration and the resulting magnetization at large distances, the dissipation leading 
to the emission can occur either in shocks or in magnetic reconnection. In both cases, 
the dominant radiative process should be synchrotron emission. Only in the first case 
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(internal shocks), detailed simulations coupling a dynamical calculation with a de¬ 
tailed radiative model are available. To have a good agreement between the observed 
spectrum and the predicted one, detailed modeling is needed, where the strength and 
structure of the magnetic field play a crucial role: moderately efficient inverse Comp¬ 
ton scatterings in the Klein-Nishina regime are needed, which requires a weak field, 
and a decay of the magnetic field far from the shock front is also probably required. 

An alternative is to explain the whole soft y-ray emission by a dissipative photo¬ 
sphere. Its nature, however, must then be elucidated. The magnetic field could again 
play a vital role, via the magnetic reconnection below the photosphere, which is a nat¬ 
ural candidate. Further progress can come from more observations over a broad spec¬ 
tral range with time-resolved spectroscopy, additional and firmer polarization mea¬ 
surements, and improvements in the modeling of the expected spectrum and spectral 
evolution in each model. Unfortunately, it remains limited by the current knowledge 
of the microphysics (structure of the magnetic field, particle acceleration) in mildly 
relativistic shocks and magnetic reconnection (Sironi et al.|2015| Kagan et al.|2015|l. 


5 Magnetic Fields in the Afterglow 

Eventually, the GRB outflow is decelerated by the external medium. It drives a strong 
relativistic blast wave - the afterglow (or external forward) shock - into the surround¬ 
ing medium. It transfers most of its energy to the shocked external medium (via pdV 
work across the contact discontinuity that separates them) at a distance R ( \ec from 
the central source - the deceleration radius. Radiation from AVc reaches the ob¬ 
server at the deceleration time, T^ ec . At R > R^ ec the original outflow composition 
no longer affects the dynamics (or emission) of the afterglow shock. However, the 
outflow magnetization can greatly affect the reverse shock (or external reverse shock, 
as it is formed due to the interaction with the external medium), whose strength and 
emission can be greatly suppressed if the outflow is highly magnetized, <7 (/Gee) > 1. 


5.1 The Afterglow Emission 

The dominant emission mechanism in the afterglow is thought to be synchrotron ra¬ 
diation, which is produced by relativistic electrons accelerated at the afterglow shock 
that gyrate in the magnetic fields within the shocked external medium^] Such a syn¬ 
chrotron origin of the afterglow emission is strongly supported by the detection of 
linear polarization at the level of ~ 1-3% in several optical or NIR afterglows (see 
§5.2| i, and by the shape of the broadband spectrum, which consists of several power 
law segments that smoothly join at some typical break frequencies. Figure [10] shows 
the possible resulting afterglow spectra. Broadband (radio to y-ray) afterglows fit 
such synchrotron spectra far better than the prompt emission. The broad and mostly 

11 In an alternative scenario, the afterglow emission is dominated at early times by the contribution of a 
long-lived reverse shock I Uhm and Beloborodov 2007 Genet et al. 2 007|, which allow s to reproduce more 
easily the observed diversity and variability, such as X-ray plateaus (Uhm et al.|2 012 Hascoet et al.|2014| 
or X-ray flares (Hascoet et al.|20l5) , though in this scenario a transition to forward shock dominance is 
expected at late times but not observed. 
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featureless smoothly broken power-law shapes of afterglow spectra evolve and fade 
more slowly over time, and have characteristic frequencies that vary as a power law 


with time, roughly according to the theoretically expected power-law indices (Sari 


et al.|1998||Granot and Sari|2002| l. Synchrotron self-Compton (SSC) - the inverse- 


Compton scattering of the synchrotron photons to (much) higher energies by the same 
population of relativistic electrons that emits the synchrotron photons - can some¬ 
times dominate the afterglow flux in the X-rays (|Sari and Esin|200~Tj |Harrison et al.| 


2001), and may affect the synchrotron emission by increasing the electron cooling. 


Relativistic collisionless shock physics (e.g., how they amplify the magnetic field 
and accelerate a non-thermal population of relativistic particles) are still not well 
understood from first principles (e.g. Sironi et al.;l 2015| l. Thus, simple assumptions 
are usually made that conveniently parameterize our ignorance. The electrons are as¬ 
sumed to be (instantly) shock-accelerated into a power-law distribution of energies, 
dNfly. rc y e p for y e > y m , and then cool both adiabatically and due to radiative 
losses| l2 | l he relativistic electrons are assumed to hold a fraction £,, of the internal 
energy immediately behind the shock, while the magnetic field is assumed to hold a 
fraction £b of the internal energy everywhere in the shocked region. Both the temporal 
and spectral indices depend on the power law index p of the electron energy distribu¬ 
tion. The temporal index (i.e. the rate of flux decay) also depends on the circumburst 
density profile, which is parameterized in Figure 10 as a power law of index k with 
the distance R from the central source, p ext R , with k = 0 and k = 2, respectively, 
corresponding to an ISM and a stellar wind - WIND. The temporal index can also be 
affected by other factors, such as energy losses or injection into the afterglow shock, 
the afterglow jet angular structure and the viewing angle relative to the jet symmetry 
axis, or time evolution of the shock microphysics parameters E e and/or t'n. 


5.2 Polarization: Afterglow and Reverse Shock 


The detection of linear polarization of a few percent in the optical and NIR afterglow 
of several GRBs (see Covino et al. |2004| and references therein) was considered as 
a confirmation that synchrotron radiation is the dominant afterglow emission mech¬ 
anism. The synchrotron emission from a fluid element with a locally uniform mag¬ 
netic field is linearly polarized in the direction perpendicular to the projection of the 
magnetic field onto the plane normal to the wave vector. Since the source moves rel¬ 
ativistic ally, one must account for aberration of light when calculating the observed 
local direction of polarization. Figure [IT] shows the predicted local polarization map 
from emission by an ultra-relativistic expanding shell, for two different magnetic field 
structures: a magnetic field that is random within the plane normal to the radial direc¬ 


tion (left panel) as could be expected from a shock-produced field (e.g., Medvedev 


and Loeb|1999]), and an ordered magnetic field normal to the radial direction (right 


panel ; as could be expected in the prompt or reverse-shock emission for a magnetic 
field coherent on angular scales 1/F that is advected from the central source). 


12 It is usually also further assumed that practically all of the electrons take part in this acceleration 
process and form such a non-thermal (power-law) distribution, leaving no thermal component (which is 
not at all clear or justified; e.g. Eichler and Waxman|2005|. 
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V in Hz 


Fig. 10 The afterglow synchrotron spectrum, calculated for the Blandford and McKee (1976) spherical 
self-similar solution, under standard assumptions, using the accurate form of the synchrotron spectral 
emissivity and integration over the emission from the whole volume of shocked material behind the for¬ 
ward (afterglow) shock (for details see |Granot and Sari|2002| from which this figure is taken). The different 
panels show the five possible broad band spectra of the afterglow synchrotron emission, each correspond¬ 
ing to a different ordering of the spectral break frequencies. Each spectrum consists of several power-law 
segments (PLSs; each shown with a different color and labeled by a different letter A-H) that smoothly 
join at the break frequencies (numbered 1-11). The broken power-law spectrum, which consists of the 
asymptotic PLSs that abruptly join at the break frequencies (and is widely used in the literature), is shown 
for comparison. Most PLSs appear in more than one of the five different broad band spectra. Indicated next 
to the arrows are the temporal scaling of the break frequencies and the flux density at the different PLSs, 
for a uniform (ISM; k — 0) and stellar wind (WIND; k = 2) external density profile. 
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Fig. 11 The predicted polarization map for synchrotron emission from a thin spherical ultra-relativistic 
shell expanding with a Lorentz factor f>l. The double-sided arrows show the direction of the linear 
polarization (the wave electric vector), while their length depends monotonically on the polarized intensity 
(in a non-trivial way, for display purposes). The circle indicates an angle oi l ff around the line of sight to 
the central source, and contains the region responsible for most of the observed flux. Left: for a magnetic 
held that is random within the plane of the shell (normal to the radial direction), for which the polarization 
direction always points at the center of the image, where the polarization vanishes (due to symmetry 
consideration). Right: for an ordered magnetic held within the plane of the shell that is coherent over 
angular scales 1/V jGranot and Konigl 2 003} . In this case the direction of the ordered magnetic held 
clearly breaks the symmetry around the center of the image, resulting in a large net polarization. For 
simplicity, the map is for a constant emission radius, rather than for a constant photon arrival time. 


The afterglow image is almost always unresolved, so we can only measure the 
(weighted) average polarization over the whole image. Therefore, a shock produced 
magnetic held that is symmetric about the shock normal will procure no net polariza¬ 
tion for a spherical how (as in this case the polarization pattern across the image is 
symmetric around its center, and the polarization averages out to zero when summed 
over the the whole image). For a shock-produced magnetic held, one thus needs to 
break this symmetry of the emission to produce net polarization. A simple and nat¬ 


ural way of doing this is considering a jet, or narrowly collimated outflow (e.g.. Sari 


1999} Ghisellini and Lazzati|l999). In this picture a jet geometry together with a line 


of sight that is not along the jet symmetry axis (but still within the jet aperture, in or¬ 
der to see the prompt GRB) is needed to break the symmetry of the afterglow image 
around our line of sight. Other models for afterglow polarization include a magnetic 
held that is coherent over patches of a size comparable to that of causally connected 
regions (Gruzinov and Waxma n|1999[ ), polarization that is induced by microlensing 
( |Loeb and Perna|1998| l or by scintillations in the radio ( |Medvedev and Loeb|1 999), 
a small ordered magnetic held component originating from the circumburst medium 
( Granot and KonigT|2 003 ), clumps in the external medium ( |Granot and Konigl 2003), 
or a very inhomogeneous jet angular structure - a patchy shell with hot spots (Granot 


and Konigl 2003; Nakar and Qren|2004l, The many possible causes of polarization. 


and the degeneracy with other factors makes it difficult to robustly determine the 
magnetic held stmcture in the emitting region from afterglow polarization measure- 
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ments. Nonetheless, a high degree of linear polarization with a stable position angle 
is hard to produce without a magnetic field that is ordered on large scales. 

The reverse shock has two main observational signatures: a sharply-peak “optical 
flash” (e.g., Akerlof et al.|1999]l on a timesc ale comparable to the prompt GRB Tgo, 
and a “radio flare” (e.g., Kulkarni et al. 1999| Frail et al. 2000||Berger et al. 2003| > that 
peaks on a timescale of a day or so after the GRB. In such cases, if the relatively bright 
observed emission is indeed from the reverse shock, this implies that the outflow was 
not strongly magnetized near the deceleration epoch, (7(/Gee) 1. Moreover, the 
polarization properties of the synchrotron emission from the reverse shock provide 
a powerful and unique probe for the magnetic field structure in the original outflow. 
Early optical polarization measurements from T > 7j ec have finally been obtained in 
the last eight years or so ( [Mundell et al.|2007a||2013]|Steele et al.|2009| l. 

On the one hand, there is a strict upper limit on the degree of linear polarization 
from GRB 060418 of P < 8% (2-cr) at T = 203 s after the GRB trigger, while the 
deceleration time suggested by the early optical lightcurve of this GRB is 7j ec = 
153 ± 10 s (Molinari et al. 2007| and its prompt emission lasted only Tgo = 52 ± 1 s. 
The fact that T^ ec /Tgo ss 3 suggest a “thin shell” in this case, which is consistent with 
a moerate magnetization (cr(Rdec) “Cl) that allows a strong reverse shock with bright 
emission. However, the polarization is fairly low near 7j ec = 153 ± 10 s (at T = 203 s), 
which suggests that either the reverse shock emission even near its peak is for some 
reason greatly sub-dominant compared to the (very weakly polarized) forward shock 
emission, or more likely that in this GRB there is hardly any ordered magnetic field 
in the ejecta on angular scales >1/1" that cover most of the visible region. 

On the other hand, GRB 090102 had a prompt duration of Tgo = 21 & and an 
optical linear polarization of P = 10.2 ± 1.3% in a 60 s exposure starting at T = 161 s 
after the trigger time (Steele et alT| |2009 1. Its optical lightcurve shows a power-law 
decay F v <^t~ a with a = 1.50±0.06 from T ~ 40 s to T ~ 1000 s and then flattens to 
a = 0.97 ±0.03 (Gendre et al. 2010). This suggests a deceleration time 7j ec < 40 s, 
well before the polarization measurement. However, the optical emission may be 
dominated by the reverse shock up to the break time of ~ 1000 s, and in particular 
during the polarization measurement. In the latter case this might possibly explain 
the measured polarization as arising from an ordered magnetic field component in the 
ejecta, though a purely ordered field on the scale of the whole emitting region (angular 
scale >1/1”) would produce a significantly larger polarization of several tens of 
percent ( |Granot and Konigl||2003| |Granot||2003 | | Lyutikov et al.||2003| |Nakar et al.| 
2003), which would suggest either a smaller magnetic field coherence length or a 


dominant contribution from the much less polarized external forward shock emission. 

Finally, GRB 120308A that lasted Tgo ~ 100 s (between T sa —30 s and T « 70 s 
post-trigger) showed an optical linear polarization of P = 28 ± 4% in an exposure 
between T = 240 s and 323 s (|Mundell et ak| ||2013[ ), which gradually decreased to 
P = 16 )j% over the next ten minutes, while keeping an approximately constant po¬ 
sition angle (to within an accuracy of about 15°). The optical lightcurve peaked at 
around T ~ 300 s, during the time bin in which the largest polarization was mea¬ 
sured, and subsequently decayed, with a possible transition from reverse to forward 
shock domination of the optical emission around ~ 1000 s. This strongly suggests 
the presence of a large-scale ordered magnetic field in the original GRB ejecta. 
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Observations of radio flares at roughly a day after the GRB have so far produced 
no detection of polarization. However, these observations have enabled to set strict 
upper limits on a possible linear or circular polarization (Granot and Taylor f2005j >. 
The strictest limits are for GRB 991216, for which the 3-ct upper limits on the linear 
and circular polarization are Pji n <7% and P c ; rc < 9%, respectively. These limits 
provide interesting constraints on existing GRB models (Granot and Taylor| j2005| >, 
and in particular are hard to reconcile with a predominantly ordered toroidal magnetic 
field in the GRB outflow together with a “structured” jet, where the energy per solid 
angle drops as the inverse square of the angle from the jet axis, as is expected in some 
models in which the outflow is Poynting flux dominated. 

Recently, the detection of circular polarization was reported in the optical after¬ 
glow of GRB 121024A by Wiersema et alT] P014| . In particular, they measured a 
circular polarization of / J c j lc = 0.61 ±0.13% at T = 0.15 days after the GRB. The 
linear polarization during that time was f\ m ~ 4% implying a circular to linear polar¬ 
ization ratio of P c i rc /P\i n ~ 0.15. A very recent detailed study that examined different 
assumptions for the magnetic field configuration, jet geometry and electron pitch- 
angle distribution (Nava et al.j2015) concluded that such a relatively high P lrc /Pm 
ratio cannot be produced by synchrotron emission from the afterglow (i.e. forward 
external) shock, which suggests an alternative origin. 


5.3 Maximum Synchrotron Photon Energy 


Since the afterglow emission from the shocked external medium is independent of 
the outflow composition, it can much more “cleanly” probe the physics of relativistic 
collisionless shocks, and serve as a testbed for how the very weak upstream magnetic 
fields in the pristine surrounding medium are amplified in the shock, and how the 
particles are accelerated in this shock and radiate in the downstream magnetic field. 

A recent challenge to the standard synchrotron afterglow scenario was raised by 
the exceptional GRB 130427A (Ackermann et al. 2014| |Perley et al. 2014| Kou- 


veliotou et al.||2013] |Maselli et al.||2014 !. This was a very energetic GRB (with an 


isotropic equivalent y-ray energy of Ey ; so = 1.4 x 10 54 erg) and it occured relatively 
nearby (at redshift z = 0.34). Therefore, it was extremely bright and in particular, 
it was detected by Fermi-LAI (100 MeV - 100 GeV) for nearly a day, including a 
95 GeV photon several minutes after the burst and a 32 GeV photon after 9 hours 
(Ackermann et al7| |2014) >. Altogether, this GRB has a large number of high-energy 
photons that clearly violate ( |Ackermann et al.|2014[ > the maximum synchrotron pho¬ 
ton energy limit] 13 1 


^syn,max 


(l±z) a 


= 3.5 


1 +z) 100 


GeV 


(18) 


which is obtained by equating the electron acceleration time to its synchrotron cool¬ 
ing time, assuming that it is accelerated and radiates its synchrotron emission in same 


13 The exact numerical coefficient depends on the exact assumptions, and in particular on whether the 
acceleration time is assumed to be a fraction of or a complete Larmor gyration time, which is in any case a 
very fast acceleration, and arguably even unrealistically so. Here a » 1/137 is the fine structure constant. 
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Fig. 12 The optical to GeV spectrum of GRB 130427A fit with the afterglow synchrotron model ofj Granot 
|and Sari| l |2002f Broadband SEDs are shown during the first ( top-panel ) and the second ( bottom-panel ) 
NuSTAR epochs. The Fermi /LAT upper-limits are shown as arrows and the extrapolation of the LAT flux 
light curve is shown as a dashed magenta cross (only during the first epoch). The second epoch ( bottom- 
pane /) is fit with a power law (black lines ); the fit to the first epoch is scaled down and superposed on the 
second epoch data for comparison (in gray). (This figure is taken from |Kouveliotou et al.|2013] ) 


same magnetic field strength. This has been argued in order to rule out an afterglow 
synchrotron origin of the late-time high-energy LAT photons, and in particular moti¬ 


vated suggestions for an origin in a distinct SSC spectral component (e.g. Ackermann 
|et al.|2014[|Tam et al.|2013l|Liu et al.|2013||Fan et al.|2013b| >. 

However, Kouveliotou et al. (20131 have shown that the optical to GeV spectrum 
is consistent with a single spectral component that very nicely matches the expecta¬ 
tion for synchrotron afterglow emission (Granot and Sa ri|2002 see Fig.[X0|). Figure[l2| 
shows their spectral fit, and demonstrates that there is hardly any room for a distinct 
(SSC) spectral component to dominate the observed flux at the highest LAT energies 
(above several GeV or so, as is needed to avoid violating T’syn.max)- This conclusion 
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is strengthened by strict upper limits on the >100 GeV flux measured by VERITAS 
at three different epochs near the first NuSTAR observation (Aliu et al .[2014j i. 

Therefore, this comprises very compelling evidence for a genuine violation of the 
£syn,max in this case, which is much harder to circumvent compared with previous 
Fermi/L AT GRBs ( |Abdo et al.||2009[ |Piran and Nakar||2010[ | Atwood et al.||2013] >. 
Thus, one should start to seriously consider how this limit can indeed be violated. A 
possible solution may lie in relaxing the assumption of a uniform magnetic field and 
instead allowing for a lower magnetic field acceleration region and a higher magnetic 
field synchrotron radiation region (e.g., |Lyutikov|20r0] |KumaT et al. 2012). Such a 
situation might arise for diffusive shock acceleration (Fermi Type I) if the tangled 
shock-amplified magnetic field decays on a short length-scale behind the shock front 
(where most of the high-energy radiation is emitted), while the highest energy elec¬ 
trons are accelerated in the lower magnetic field further downstream (Kumar et al. 


2012|). In fact, such a scenario has recently been suggested (Lemoine|2013 Lemoine 


et al. 20131, and also shown to significantly alleviate the previously very high y-ray 


radiative efficiencies inferred for most Swift GRBs (Beniamini et al. 2015). 


6 Conclusions 


In this review we have demonstrated that magnetic fields clearly play a vital role in 
GRBs, practically in every important aspect of this phenomenon. Here we briefly 
summarize some of our main conclusions, and stress both what was learned so far, as 
well as what still needs to be carefully studied. Following the main text, the discussion 
proceeds more or less in order of increasing distance from the central source. 

Magnetic fields most likely play a crutial role in the launching of GRB jets. More¬ 
over, hydromagnetic jet launching implies dynamically strong magnetic fields near 
the central source, which can naturally help avoid an excessive baryonic loading into 
the jet. Such a low baryon loading is essential in order for the jet to be able to reach 
sufficiently large Lorentz factors (F > 100 — 500) that are inferred from prompt GRB 
observations. The jet launching definitely requires many further detailed studies, also 
(semi-) analytic, but mainly numerical studies. The latter are, however, involved as 
they require general-relativistic magneto-hydrodynamics (GRMHD) codes, coupled 
with neutrino, plasma and radiation physics. Therefore, is it likely to take many years 
before such studies will provide definitive answers. 

Millisecond-magnetar models for the GRB central engine, discussed in jj2] face 
serious challenges that still must be overcome. Models where the GRB arises from the 
delayed collapse of a supramassive millisecond magnetar (such as the “time reversal 
model”; Rezzolla and Kumar|2015[|Ciolfi an d Siegel 2015) face serious difficulty as 
this requires the formation of a disk during this collapse, which was recently argued 
to not be possible ( [Margalit et al.|20l5] l. Producing both the prompt GRB emission 
and the X-ray plateau observed by Swift in long GRBs is challenging for millisecond- 
magnetar models. Models invoking millisecond magnetars in short GRBs also face 
many problems, such as how to produce the short GRB, or hide the huge amount of 
of rotational energy (a few x 10 52 erg) that is injected into the afterglow shock, while 
short GRB afterglows are very dim. More detailed studies of the relevant physics, 
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both analytic and numerical, are definitely needed in order to produce more robust 
and realistic predictions that could be tested in more detail against the relevant obser¬ 
vations. A relatively simple example is that a highly-relativistic pulsar-like magnetar 
MHD wind is invoked to explain the X-ray plateaus observed by Swift. However, 
most studies consider only the forward shock emission and ignore the emission from 
the MHD wind itself (see, however, Dai|2004) , which may involve magnetic recon¬ 
nection and the resulting particle acceleration at the termination shock of the pulsar 
wind, similar to pulsar wind nebulae (see, e.g., Lyubarsky and Kirk|2001| Lyubarsky 


2003 Sironi and Spitkovsky 2014 ). The birth process of a millisecond magnetar, ei¬ 


ther in the core-collapse of a massive star or in a binary neutron star merger, and 
how it evolves in these messy environments, is obviously very involved but there is 
definitely a lot of room for improvement and new studies on such systems. 

The GRB outflow composition, and in particular its degree of magnetization is an 
important open question. On the one hand, in (3.1 we have shown very compelling 
evidence based on the theoretical jet propagation time in the progenitor star of a long 
GRB and the observed GRB duration distribution, that the jet appears to have modest 
magnetization (a < 1) throughout most of the time it takes for it to bore its way out of 
the star. Moreover, in ; |4.5| we argued based on the prompt GRB observations (both y- 
ray detections and upper limits in optical and soft X-ray) that the emission region has 
a low magnetization. On the other hand, hydromagnetic jet launching is much more 
promising than a pure thermal one, and can naturally help avoid excessive baryon 
loading, suggesting a high magnetization near the source (Oq 1). A similar con¬ 
clusion is also strongly suggested from the fact that quasi-thermal components in the 
prompt GRB spectrum are typically highly sub-dominant (see (j4]). Taking these two 
lines of evidence together, it appears that the magnetization significantly decreases 
as the jet propagates from the source out to large distances. This can occur through 
two main channels (or some combination of the two): (i) very efficient conversion 
of magnetic to kinetic energy that leaves a low magnetization in the emission region 
and allows for efficient internal shocks, through first quasi-steady and then impulsive 
magnetic acceleration (see Fig. [6] [Granot et al.|201 1 Granot 2012b), or (ii) strong 
magnetic reconnection that converts magnetic energy into thermal energy and accel¬ 
erates particles, and yet somehow leaves a sufficiently low magnetization where these 
particles radiate most of their energy. 

Magnetic fields can also play a very important role in the energy dissipation that 
leads to the GRB emission. On the one hand, sufficiently low magnetization (a < 1 
or even a < 0.1) is needed for efficient energy dissipation in internal shocks. On 
the other hand, while a large magnetization (<7 > 1) effectively suppresses internal 
shocks, or the reverse shock, it can lead to significant dissipation through magnetic 
reconnection under appropriate conditions (e.g. if the source ejects outflow with a 
magnetic field that occasionally changes its polarity, or through certain instabilities). 

Another vital role of magnetic fields in GRBs is in particle acceleration. Within 
the outflow the magnetization can be high, allowing efficient magnetic reconnection, 
which can directly convert a good fraction of the dissipated magnetic energy into the 
random motions of the particles that it accelerates in this process. The exact magnetic 
field strength and geometry throughout the emission region also greatly affect the 
radiation of the accelerated particles, and therefore their observable signatures. In 
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particular, magnetic fields are vital to the synchrotron emission, which dominates in 
the afterglow and also likely plays a key role in the prompt GRB emission. 

Both optically thin internal shocks as well as the external forward and reverse 
shocks are collisionless, and mediated by electromagnetic fields (through collective 
plasma interactions). Moreover, while the outflow itself typically has a large enough 
magnetic field for efficient synchrotron radiation, the external medium has a very 
low magnetization that must be significantly amplified at the afterglow shock front in 
order to produce the observed afterglow emission. Therefore, magnetic field amplifi¬ 
cation in relativistic collisionless shocks and its possible decay behind the shock are 
vital for understanding the afterglow physics and interpreting afterglow observations. 
Moreover they can strongly affect the particle acceleration in the afterglow shock, 
and the resulting afterglow emission. Such physics might hold the key to unravel a 
puzzle (see (5.3 i arising from observations of the very bright and relatively nearby 
long GRB 130427A, which show an apparent violation of the maximum synchrotron 
photon energy limit, T S yn.max- 

The different roles of magnetic fields in GRBs are numerous and diverse. Many of 
them are only starting to be understood, while new roles are still being occasionally 
discovered. Since magnetic fields appear in almost all aspects of GRB physics, future 
studies of their properties and effects are likely to greatly improve our understanding 
of GRBs, and lead to fundamental progress. 
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